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Abstract. Now that the low temperature properties of quantum-mech- 
anical many-body systems (bosons) at low density, p, can be examined 
experimentally it is appropriate to revisit some of the formulas deduced 
by many authors 4-5 decades ago. For systems with repulsive (i.e. posi- 
tive) interaction potentials the experimental low temperature state and 
the ground state are effectively synonymous - and this fact is used in 
all modeling. In such cases, the leading term in the energy/particle is 
2-Kh^ap/m where a is the scattering length of the two-body potential. 
Owing to the delicate and peculiar nature of bosonic correlations (such 
as the strange N"^^^ law for charged bosons), four decades of research 
failed to establish this plausible formula rigorously. The only previous 
lower bound for the energy was found by Dyson in 1957, but it was 
14 times too small. The correct asymptotic formula has recently been 
obtained by us and this work will be presented. The reason behind 
the mathematical difhculties will be emphasized. A different formula, 
postulated as late as 1971 by Schick, holds in two-dimensions and this, 
too, will be shown to be correct. With the aid of the methodology de- 
veloped to prove the lower bound for the homogeneous gas, two other 
problems have been successfully addressed. One is the proof by us that 
the Gross-Pitaevskii equation correctly describes the ground state in the 
'traps' actually used in the experiments. For this system it is also pos- 
sible to prove complete Bose condensation, as we have shown. Another 
topic is a proof that Foldy's 1961 theory of a high density Bose gas of 
charged particles correctly describes its ground state energy. All of this 
is quite recent work and it is hoped that the mathematical methodology 
might be useful, ultimately, to solve more complex problems connected 
with these interesting systems. 



Foreword 

At the conference "Contemporary Developments in Mathematics" , hosted 
by the MIT and Harvard University Mathematics Departments, November 
16-17, 2001, one of us (E.H.L.) contributed a talk with the title "The Bose 
gas: A subtle many-body problem". This talk covered material by all the 
authors listed above. This contribution is a much expanded version of the 
talk and of [L3j . 
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1. Introduction 

Schrodinger's equation of 1926 defined a new mechanics whose Hamilton- 
ian is based on classical mechanics, but whose consequences are sometimes 
non-intuitive from the classical point of view. One of the most extreme cases 
is the behavior of the ground (= lowest energy) state of a many-body system 
of particles. Since the ground state function "^{xi, xn) is automatically 
symmetric in the coordinates Xj £ of the N particles, we are dealing 
necessarily with 'bosons '. If we imposed the Pauli exclusion principle (anti- 
symmetry) instead, appropriate for electrons, the outcome would look much 
more natural and, oddly, more classical. Indeed, the Pauli principle is es- 
sential for understanding the stability of the ordinary matter that surrounds 
us. 

Recent experiments have confirmed some of the bizarre properties of 
bosons close to their ground state, but the theoretical ideas go back to 
the 1940's - 1960's. The first sophisticated analysis of a gas or liquid of 
interacting bosons is due to Bogolubov in 1947. His approximate theory as 
amplified by others, is supposed to be exact in certain limiting cases, and 
some of those cases have now been verified rigorously (for the ground state 
energy) — 3 or 4 decades after they were proposed. 

The discussion will center around four main topics. 

(1) The dilute, homogeneous Bose gas with repulsive interaction (2D 
and 3D). 

(2) Repulsive bosons in a trap (as used in recent experiments) and the 
'Gross-Pitaevskii' equation. 

(3) Bose-Einstein condensation for dilute trapped gases. 
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(4) Foldy's 'jellium' model of charged particles in a neutralizing back- 
ground. 

Note that for potentials that tend to zero at infinity 'repulsive' and 'pos- 
itive' are synonymous — in the quantum mechanical literature at least. In 
classical mechanics, in contrast, a potential that is positive but not mono- 
tonically decreasing is not called repulsive. 

The discussion below of topic 1 is based on |LYlj and |LY2j , and of topic 2 
on ILSeYl ' and ' LSeY2j . See also [LY3llLSeY3lF^ILSeY4| . The discussion 
of topic 3 is mainly taken from |LSej . but for transparency we also include 
here a section on the special case when the trap is a rectangular box. This 
case already contains the salient points, but avoids several complications 
due the the inhomogeneity of the gas in a general trap. The discussion of 
topic 4 is based on |LSoj . 

Topic 1 (3-dimensions) was the starting point and contains essential ideas. 
It is explained here in some detail and is taken, with minor modifications 
(and corrections), from |LY3j . In terms of technical complexity, however, 
the fourth topic is the most involved and can not be treated here in full 
detail. 

The interaction potential between pairs of particles in the Jellium model 
in topic 4 is the repulsive, long-range Coulomb potential, while in topics 
1-3 it is assumed to be repulsive and short range. For alkali atoms in the 
recent experiments on Bose Einstein condensation the interaction potential 
has a repulsive hard core, but also a quite deep attractive contribution of 
van der Waals type and there are many two body bound states |PSj . The 
Bose condensate seen in the experiments is thus not the true ground state 
(which would be a solid) but a metastable state. Nevertheless, it is usual to 
model this metastable state as the ground state of a system with a repulsive 
two body potential having the same scattering length as the true potential, 
and this is what we shall do. In this paper all potentials will be positive. 

2. The Dilute Bose Gas in 3D 

We consider the Hamiltonian for N bosons of mass m enclosed in a cubic 
box A of side length L and interacting by a spherically symmetric pair 
potential v{\xi — Xj\): 

N 

i?Ar = -/i^Ai+ v{\xi-Xj\). (2.1) 

1=1 l<i<j<N 

Here Xi G M^, i = 1, . . . ,N are the positions of the particles, Aj the Lapla- 
cian with respect to Xj, and we have denoted fi^ /2m by /i for short. (By 
choosing suitable units could, of course, be eliminated, but we want to 
keep track of the dependence of the energy on Planck's constant and the 
mass.) The Hamiltonian (|2.1j) operates on symmetric wave functions in 
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L'^{A^ ,dxi ■ ■ ■ dx]\f) as is appropriate for bosons. The interaction poten- 
tial will be assumed to be nonnegative and to decrease faster than at 
infinity. 

We are interested in the ground state energy Eq[N,L) of (|2.1j) in the 
thermodynamic limit when N and L tend to infinity with the density p = 
N/L^ fixed. The energy per particle in this limit is 

eo(p)= lim i?o(pi',i)/(pi'). (2.2) 

L^oo 

Our results about eo(/3) are based on estimates on Eo{N, L) for finite and 
L, which are important, e.g., for the considerations of inhomogeneous sys- 
tems in |LSeYlj . To define Eq{N, L) precisely one must specify the boundary 
conditions. These should not matter for the thermodynamic limit. To be 
on the safe side we use Neumann boundary conditions for the lower bound, 
and Dirichlet boundary conditions for the upper bound since these lead, 
respectively, to the lowest and the highest energies. 

For experiments with dilute gases the low density asymptotics of eo{p) is 
of importance. Low density means here that the mean interparticle distance, 
p-i/3 jg j^ucJi larger than the scattering length a of the potential, which is 
defined as follows. The zero energy scattering Schrodinger equation 

-2pAi; + v{r)i{j = (2.3) 

has a solution of the form, asymptotically as \x\ = r — > oo (or for all r > Rq 
if v{r) = for r > Rq), 

i^Q{x) = 1 - a/\x\ (2.4) 

(The factor 2 in (|2.3|1 comes from the reduced mass of the two particle 
problem.) Writing i^oix) = -uod^^D/l^^l this is the same as 

a = lim r - (2.5) 

r^oo UQyr) 

where uq solves the zero energy (radial) scattering equation, 

-2/iUo(r) + v{r)uQ{r) = (2.6) 

with no(0) = 0. 

An important special case is the hard core potential v{r) = oo \i r < a 
and v{r) =0 otherwise. Then the scattering length a and the radius a are 
the same. 

Our main result is a rigorous proof of the formula 

eo(/o) K. AiTfipa (2.7) 

for pa^ <^ 1, more precisely of 

Theorem 2.1 (Low density limit of the ground state energy). 

lim ^ = L (2.8) 
pa3->o Airppa 
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This formula is independent of the boundary conditions used for the defi- 
nition of eo(p) • It holds for every positive radially symmetric pair potential 
such that v{r)r'^dr < oo for some R, which guarantees a finite scattering 
length, cf. Appendix A in |LY2j . 

The genesis of an understanding of eo{p) was the pioneering work jBoj of 
Bogolubov, and in the 50's and early 60's several derivations of (|2.8|) were 
presented |HYj . |Llj . even including higher order terms: 

""^'^ - 1 + ^ipa'?/' + 8(f-V3) {pa') log{pa') + 0{pa') (2.9) 



Anppa \^ 3 

These early developments are reviewed in |L2j . They all rely on some special 
assumptions about the ground state that have never been proved, or on the 
selection of special terms from a perturbation series which likely diverges. 
The only rigorous estimates of this period were established by Dyson, who 
derived the following bounds in 1957 for a gas of hard spheres |Dlj : 



1 eo(p) l + 2yV3 

^ < (2.10) 



WV2 ~ 4.T^ppa ~ (1 - y 1/3)2 

with Y = A-npa?/?,. While the upper bound has the asymptotically correct 
form, the lower bound is off the mark by a factor of about 1/14. But for 
about 40 years this was the best lower bound available! 

Under the assumption that (|2.8|) is a correct asymptotic formula for the 
energy, we see at once that understanding it physically, much less proving it, 
is not a simple matter. Initially, the problem presents us with two lengths, 
a <C /0~i/3 at low density. However, 1)2. 8|) presents us with another length 
generated by the solution to the problem. This length is the de Broglie wave- 
length, or 'uncertainty principle' length (sometimes called 'healing length') 

4~(H-V2. (2.11) 

The reason for saying that is the de Broglie wavelength is that in the 
hard core case all the energy is kinetic (the hard core just imposes a = 
boundary condition whenever the distance between two particles is less than 
a). By the uncertainty principle, the kinetic energy is proportional to an 
inverse length squared, namely Ic- We then have the relation (since pa? is 
small) 

a < p-i/3 < 4 (2.12) 

which implies, physically, that it is impossible to localize the particles rela- 
tive to each other (even though p is small). Bosons in their ground state are 
therefore 'smeared out' over distances large compared to the mean particle 
distance and their individuality is entirely lost. They cannot be localized 
with respect to each other without changing the kinetic energy enormously. 

Fermions, on the other hand, prefer to sit in 'private rooms', i.e., 4 is 
never bigger than p~^^' by a fixed factor. In this respect the quantum 
nature of bosons is much more pronounced than for fermions. 
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Since (|2.8j) is a basic result about the Bose gas it is clearly important 
to derive it rigorously and in reasonable generality, in particular for more 
general cases than hard spheres. The question immediately arises for which 
interaction potentials one may expect it to be true. A notable fact is that 
it not true for all v with a > 0, since there are two body potentials with 
positive scattering length that allow many body bound states. (There are 
even such potentials without two body bound states but with three body 
bound states tBaj .) For such potentials (|2.8j) is clearly false. Our proof, 
presented in the sequel, works for nonnegative v, but we conjecture that 
(|2.8() holds if a > and v has no A^-body bound states for any A^. The 
lower bound is, of course, the hardest part, but the upper bound is not 
altogether trivial either. 

Before we start with the estimates a simple computation and some heuris- 
tics may be helpful to make 1)2. 8|) plausible and motivate the formal proofs. 

With ipQ the zero energy scattering solution, partial integration, using 
(tO) and ((231), gives, for R > Rq, 

/ {2^|VV'o|^ + v\'4^o\'^}dx = STTfia (l — j Sivfia for R ^ oo. 

J\x\<R ^ 

(2.13) 

Moreover, for positive interaction potentials the scattering solution mini- 
mizes the quadratic form in 1)2. 13(1 for each R> Rq with the boundary con- 
dition V'o(|2^| = R) = (1 — a/i?). Hence the energy Eq(2, L) of two particles 
in a large box, i.e., L ^ a, is approximately Sir fj,a / L"^ . If the gas is suffi- 
ciently dilute it is not unreasonable to expect that the energy is essentially 
a sum of all such two particle contributions. Since there are A^(A^ — l)/2 
pairs, we are thus lead to Eo{N, L) ^ 4TTnaN{N - 1)/L^, which gives 1)2. 8(1 
in the thermodynamic limit. 

This simple heuristics is far from a rigorous proof, however, especially for 
the lower bound. In fact, it is rather remarkable that the same asymptotic 
formula holds both for 'soft' interaction potentials, where perturbation the- 
ory can be expected to be a good approximation, and potentials like hard 
spheres where this is not so. In the former case the ground state is approxi- 
mately the constant function and the energy is mostly potential: According 
to perturbation theory £'o(A^, L) ~ A^(A^ — 1)/(2L^) J v{\x\)dx. In particular 
it is independent of /i, i.e. of Planck's constant and mass. Since, however, 
/ v{\x\)dx is the first Born approximation to 87r ^ua (note that a depends on 
/i!), this is not in conflict with (|2.8|) . For 'hard' potentials on the other hand, 
the ground state is highly correlated, i.e., it is far from being a product of 
single particle states. The energy is here mostly kinetic, because the wave 
function is very small where the potential is large. These two quite different 
regimes, the potential energy dominated one and the kinetic energy dom- 
inated one, cannot be distinguished by the low density asymptotics of the 
energy. Whether they behave differently with respect to other phenomena, 
e.g., Bose-Einstein condensation, is not known at present. 
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Bogolubov's analysis |Boj presupposes the existence of Bose-Einstein con- 
densation. Nevertheless, it is correct (for the energy) for the one-dimensional 
delta-function Bose gas |LLj . despite the fact that there is (presumably) no 
condensation in that case |PiStj . It turns out that BE condensation is not 
really needed in order to understand the energy. As we shall see, 'global' 
condensation can be replaced by a 'local' condensation on boxes whose size 
is independent of L. It is this crucial understanding that enables us to prove 
Theorem 1.1 without having to decide about BE condensation. 

An important idea of Dyson was to transform the hard sphere potential 
into a soft potential at the cost of sacrificing the kinetic energy, i.e., effec- 
tively to move from one regime to the other. We shall make use of this 
idea in our proof of the lower bound below. But first we discuss the simpler 
upper bound, which relies on other ideas from Dyson's beautiful paper |Dlj . 



2.1. Upper Bound. The following generalization of Dyson's upper bound 
holds jL^eYlj . f^ : 



Theorem 2.2 (Upper bound). Define pi = {N-l)/L^ andb = (47rpi/3)"i/3. 
For nonnegative potentials v and b > a the ground state energy of i2.1\) with 
periodic boundary conditions satisfies 

Eo{N, L)/N < "y^^V • (2-14) 

V- ~ b) 

For Dirichlet boundary conditions the estimate holds with (const. )/L^ added 
to the right side. Thus in the thermodynamic limit and for all boundary 
conditions 

e„(p) ^ i-rv^ + rV3-ir ^^^^^^ 



47r/ipa - (1 - yi/3)8 

provided Y = 'iir < 1. 

Remark. The bound ()2.14|) holds for potentials with infinite range, provided 
b > a. For potentials of finite range Rq it can be improved for b > Rq to 

1 _ (il)'^ + 1 (a)^ 
Eo{N, L)/N < Anpp.a ^-^ (2.16) 

(1-f) 

Proof. We first remark that the expectation value of 1)2. 1() with any trial 
wave function gives an upper bound to the bosonic ground state energy, 
even if the trial function is not symmetric under permutations of the vari- 
ables. The reason is that an absolute ground state of the elliptic differential 
operator ()2.1|) (i.e., a ground state without symmetry requirement) is a non- 
negative function which can be be symmetrized without changing the energy 
because 1)2. 1|) is symmetric under permutations. In other words, the absolute 
ground state energy is the same as the bosonic ground state energy. 

Following |Dlj we choose a trial function of the following form 

^'(fi, . . .,xn) = Fi{xi) ■ F2{xi,X2) ■ • •FAr(fi, . . .,xn)- (2.17) 
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More specifically, Fi = 1 and Fi depends only on the distance of Xi to 
its nearest neighbor among the the points xi, . . . ,Xi-i (taking the periodic 
boundary into account): 

Fi{xi, . . . ,Xi) = f{ti), U = m.m{\xi - Xj\,j = 1, . . . ,i - 1) , (2.18) 

with a function / satisfying 

< / < 1, /' > 0. (2.19) 

The intuition behind the ansatz (|2.17|) is that the particles are inserted into 
the system one at the time, taking into account the particles previously in- 
serted. While such a wave function cannot reproduce all correlations present 
in the true ground state, it turns out to capture the leading term in the en- 
ergy for dilute gases. The form ()2.17|1 is computationally easier to handle 
than an ansatz of the type Y\i^jf{\xi — Xj\), which might appear more 
natural in view of the heuristic remarks after Eq. ()2.13() . 
The function / is chosen to be 

;,,) f AM/AW forO<r<6. 
1 for r > b, 



with foir) = uo{r)/r the zero energy scattering solution defined by (|2.6j) . 
The estimates (|2.14|l and (|2.16|) are obtained by somewhat lengthy compu- 



tations similar as in |Dlj . but making use of (|2.13|) . For details we refer to 
|T.SeY1j and [Sel] . 

A test wave function with Dirichlet boundary condition may be obtained 
by localizing the wave function 1)2. 17() on the length scale L. The energy 
cost per particle for this is (const. )/L^. ■ 

2.2. Lower Bound. It was explained previously in this section why the 
lower bound for the bosonic ground state energy of 1)2. 1|) is not easy to 
obtain. The three different length scales 1)2.12(1 for bosons will play a role in 
the proof below. 

• The scattering length a. 

• The mean particle distance p~^^^. 

• The 'uncertainty principle length' £c, defined by fii~'^ = eo{p), i.e., 
4 ~ (pa)-V2. 

Our lower bound for eo(/o) is as follows. 

Theorem 2.3 (Lower bound in the thermodynamic limit). For a 

positive potential v with finite range and Y small enough 

^ol^>(l_Cyi/i7) (2.21) 
4:7rppa 

with C a constant. If v does not have finite range, but decreases faster than 
1/r^ (more precisely, v{r)r'^dr < oo for some R) then an analogous 
bound to \2.21\) holds, but with CY^^I^'^ replaced by o(l) asY^O. 
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It should be noted right away that the error term —CY^I^'^ in (|2.2H) is 
of no fundamental significance and is not believed to reflect the true state 
of affairs. Presumably, it does not even have the right sign. We mention in 
passing that C can be taken to be 8.9 |Selj . 

As mentioned at the beginning of this section after Eq. (|2.2|) . a lower 
bound on Eq{N,L) for finite N and L is of importance for applications to 
inhomogeneous gases, and in fact we derive (|2.2H) from such a bound. We 
state it in the following way: 

Theorem 2.4 (Lower bound in a finite box). For a positive potential v 
with finite range there is a 6 > such that the ground state energy of ^2. 1]) 
with Neumann boundary conditions satisfies 

Eo{N, L) /N > A-Kiipa (l-C F^/^^) (2.22) 



for all N and L with Y < 6 and L/a > C'Y-^I^'^ . Here C and C are 
positive constants, independent of N and L. (Note that the condition on L/a 
requires in particular that N must be large enough, N > {const. )Y ^ .) 
As in Theorem \2.!^ such a bound, but possibly with a different error term 
holds also for potentials v of infinite range that decrease sufficiently fast at 
infinity. 

The first step in the proof of Theorem 12.41 is a generalization of a lemma 
of Dyson, which allows us to replace f; by a 'soft' potential, at the cost of 
sacrificing kinetic energy and increasing the effective range. 

Lemma 2.5. Let v{r) > with finite range Rq. Let U{r) > be any 
function satisfying f U{r)r'^dr < 1 and U{r) = for r < Rq. Let B C M.^ 
be star shaped with respect to (e.g. convex with £ B). Then for all 
differentiable functions tp 

[ [^IWI' + it^lVl'] >l^a [ C/IV'P. (2.23) 

Proof. Actually, 1)2. 23(1 holds with ii\Vip{x)\'^ replaced by the (smaller) radial 
kinetic energy, ^\dijj{x) / dr\'^ , and it suffices to prove the analog of p.23j) for 
the integral along each radial line with fixed angular variables. Along such 
a line we write 'iIj{x) = u{r)/r with ii(0) = 0. We consider first the special 
case when when [/ is a delta-function at some radius R> Rq, i.e., 

U{r) = ^6ir-R). (2.24) 

For such U the analog of 1)2. 23() along the radial line is 

- (u(ryr)f + i„(0N.)i^)* > I ^'^^^ 

(2.25) 

where Ri is the length of the radial line segment in B. The case Ri < R is 
trivial, because ii\dijj / dr]"^ + ^t'lV'P > 0. (Note that positivity of v is used 
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here.) If i? < i?i we consider the integral on the the left side of (|2.25j) from 
to i? instead of Ri and minimize it under the boundary condition that 
u{0) = and u{R) is a fixed constant. Since everything is homogeneous in 
u we may normahze this value to u{R) = R — a. This minimization problem 
leads to the zero energy scattering equation p.6|) . Since v is positive, the 
solution is a true minimum and not just a stationary point. 

Because v{r) =0 for r > Rq the solution, uq, satisfies UQ(r) = r — a for 
r > Rq. By partial integration, 

R 

{^l[u'f^{r)-{uo{r)/r)f + ^v{r)\uo{r)]^}dr = fia\R-a\/R > fia\R-a\'^ /R^ . 

(2.26) 

But — ap/i?^ is precisely the right side of (|2.25|) if u satisfies the normal- 
ization condition. 

This derivation of (|2.23|1 for the special case (|2.24|) implies the general 
case, because every U can be written as a superposition of (^-functions, 
U{r) = J R-'^5{r-R) U{R)R^dR, and / U{R)R^dR < 1 by assumption. ■ 

By dividing A for given points xi, . . . ,xn into Voronoi cells Bi that contain 
all points closer to Xi than to Xj with j ^ i (these cells are star shaped w.r.t. 
Xi , indeed convex) , the following corollary of Lemma 12.51 can be derived in 
the same way as the corresponding Eq. (28) in |Dlj . 

Corollary 2.6. For any U as in Lemma \2. 51 

Hn > IJ-aW (2.27) 

with 

N 

iy(fi,...,f7v) =X]^(**)' (2-28) 
1=1 

where ti is the distance of Xi to its nearest neighbor among the other points 
Xj, j = 1,...,N, i.e., 

ti{xi, . . . , xj\f) = min \xi — Xj\. (2.29) 

(Note that ti has here a slightly different meaning than in 1)2. 18() . where it 
denoted the distance to the nearest neighbor among the xj with j < i — 1.) 

Dyson considers in |Dlj a one parameter family of C/'s that is essentially 
the same as the following choice, which is convenient for the present purpose: 

I otherwise. 

We denote the corresponding interaction 1)2. 28() by Wr. For the hard core 
gas one obtains 

^(iV,L)>sup inf iuaWR{xi,...,XN) (2.31) 

R (xi,...,Xn) 
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where the infimum is over {xi, . . . ,X]\f) S with \xi — Xj\ > Rq = a, 
because of the hard core. At fixed R simple geometry gives 

with certain constants A and B. An evaluation of these constants gives 
Dyson's bound 

EiN,L)/N>,^4nfipa. (2.33) 

The main reason this method does not give a better bound is that R 
must be chosen quite big, namely of the order of the mean particle distance 
p"^/'^, in order to guarantee that the spheres of radius R around the 
points overlap. Otherwise the infimum of Wn will be zero. But large R 
means that Wr is small. It should also be noted that this method does not 
work for potentials other than hard spheres: If | allowed to be less 

than Rq, then the right side of (|2.31|) is zero because U{r) = for r < Rq. 

For these reasons we take another route. We still use Lemma 12.51 to get 
into the soft potential regime, but we do not sacrifice all the kinetic energy 
as in (|2.27l) . Instead we write, for e > 

Hn = eHn + (1 - e)HN > eTn + (1 - e)HN (2.34) 

with Tjv = — Ai and use (|2.27j) only for the part (1 — e)HN- This gives 

Hn > sTn + (1 - e)tiaWR. (2.35) 

We consider the operator on the right side from the viewpoint of first order 
perturbation theory, with eTjy as the unperturbed part, denoted Hq. 

The ground state of Hq in a box of side length L is ^o{xi, . . . ,xn) = 
L-3iV/2 and we denote expectation values in this state by (•)o. A computa- 
tion, cf. Eq. (21) in |LYlj . gives 

4vrp(l-i) > {Wr)q/N 

> Attp (1 - ^) (1 - f )^ (1 + 47rp(l - i)(i?3 _ Rl)/3)y' . 

(2.36) 

The rationale behind the various factors is as follows: (1 — ■^) comes from 
the fact that the number of pairs is A^(A^ - l)/2 and not (1 - 2R/Lf 

takes into account the fact that the particles do not interact beyond the 
boundary of A, and the last factor measures the probability to find another 
particle within the interaction range of the potential Ur for a given particle. 

The estimates (|2.36|) on the first order term look at first sight quite promis- 
ing, for if we let L — > oo, N ^ oo with p = N/ fixed, and subsequently 
take R — > oo, then (W/j)o/A^ converges to 47rp, which is just what is desired. 
But the first order result (|2.36|) is not a rigorous bound on Eq{N,L), we 
need error estimates, and these will depend on e, R and L. 
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We now recall Temple 's inequality |Tj for the expectation values of an op- 
erator H = Hq + V in the ground state (•)o of Hq. It is a simple consequence 
of the operator inequality 

{H - Eo){H - El) > (2.37) 
for the two lowest eigenvalues, Eq < Ei, of H and reads 

^0 > (^)o - ^^^I^ (2-38) 

provided Ei — {H)q > 0. Furthermore, if V > we may use Ei > Ef^^ = 

second lowest eigenvalue of Hq and replace Ei in (|2.38|) by E^^ . 
From 1)2. 36() and 1)2. 38|) we get the estimate 

Eo{N, L) 
N 

with 



>iTrnap{l-£{p,L,R,e)) (2.39) 



Ma((W^l)o - {Wr)1) 



V {WR)o{Ef^ - fia{WR)o) J ' 
To evaluate this further one may use the estimates (|2.36() and the bound 

{W'r)o<3j^^^JWr)o (2.41) 

which follows from = 3(i?^ — Rq)^^Ur together with the Cauchy-Schwarz 
inequality. A glance at the form of the error term reveals, however, that it 
is not possible here to take the thermodynamic limit L — > oo with p fixed: 
We have E^^ = etrp/L? (this is the kinetic energy of a single particle in 

the first excited state in the box), and the factor Ef'^ — pa{WR)o in the de- 
nominator in ()2.4U|) is, up to unimportant constants and lower order terms, 
~ (sL~'^ — ap^L^). Hence the denominator eventually becomes negative and 
Temple's inequality looses its validity if L is large enough. 

As a way out of this dilemma we divide the big box A into cubic cells of 
side length £ that is kept fixed as L — > cxd. The number of cells, /£^, on the 
other hand, increases with L. The N particles are distributed among these 
cells, and we use (|2.4()|) . with L replaced by i, N by the particle number, n, 
in a cell and p by n/i^, to estimate the energy in each cell with Neumann 
conditions on the boundary. For each distribution of the particles we add 
the contributions from the cells, neglecting interactions across boundaries. 
Since f > by assumption, this can only lower the energy. Finally, we 
minimize over all possible choices of the particle numbers for the various 
cells adding up to A^. The energy obtained in this way is a lower bound to 
Eq(N,L), because we are effectively allowing discontinuous test functions 
for the quadratic form given by H]\f. 



(2.40) 
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In mathematical terms, the ceh method leads to 

Eo{N,L)/N > {pf)-'miY,CnEo{n,e) (2.42) 

n>0 

where the infimum is over all choices of coefficients c„ > (relative number 
of cells containing exactly n particles), satisfying the constraints 

^c„ = l, ^Cnn = pi\ (2.43) 

n>0 n>0 

The minimization problem for the distributions of the particles among 
the cells would be easy if we knew that the ground state energy Eo(n,£) (or 
a good lower bound to it) were convex in n. Then we could immediately 
conclude that it is best to have the particles as evenly distributed among 
the boxes as possible, i.e., c„ would be zero except for the n equal to the 
integer closest to p£^. This would give 

> 47Tpap (1 - £{p, £, R, e)) (2.44) 

i.e., replacement of L in (|2.39|) by i, which is independent of L. The blow 
up of £^ for L ^ oo would thus be avoided. 

Since convexity of EQ{n,£) is not known (except in the thermodynamic 
limit) we must resort to other means to show that n = 0{pi^) in all boxes. 
The rescue comes from superadditivity of EQ{n,£), i.e., the property 

Eo{n + n',i)>Eo{n,i) + Eo{n,l) (2.45) 

which follows immediately from f > by dropping the interactions between 
the n particles and the n' particles. The bound 1)2. 45(1 implies in particular 
that for any n,p gN with n> p 

Tl 

E(n, e) > [n/p] E(p, £) > —E(p, () (2.46) 

2p 

since the largest integer [n/p] smaller than n/p is in any case > n/{2p). 

The way (|2.46|) is used is as follows: Replacing L by ^, by n and p by 
n/l^ in (|2.39)) we have for fixed R and e 

Eo{n,e) > ^n{n-l)K{n,£) (2.47) 

with a certain function K{n,i) determined by ()2.4U|) . We shall see that K 
is monotonously decreasing in n, so that if p E N and n < p then 

Eoin, I) > ^n(n - \)K{j,, I). (2.48) 

We now split the sum in (|2.42|) into two parts. For n < p we use 1)2. 48() . and 
for n > p we use (|2.46|) together with (|2.48j) for n = p. The task is thus to 
minimize 

Cnn{n + ^""^^^ - 1) (2.49) 
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subject to the constraints (|2.4H|) . Putting 

k := and t := ^ c^n < A: (2.50) 

n<p 

we have X]„>p c„,n = k — t, and since n(n — 1) is convex in n and vanishes 
for n = 0, and J2n<p^'n- — ^' expression 1)2. 49(1 is 

>t{t-l) + l{k-t){p-l). (2.51) 

We have to minimize this for 1 < t < k. Ifp>4A; the minimum is taken at 
t = k and is equal to k{k — 1). Altogether we have thus shown that 

Eo{N, L) ^ ^^^^^ K{4pf,£). (2.52) 



What remains is to take a closer look at K{Ap£^,£), which depends on 
the parameters e and R besides i, and choose the parameters in an optimal 
way. From (|2.40j) and (|2.41j) we obtain 

K{n,i) = {l-e){l-^-^f[l + fp[l-l){R^-Rl)))-' 

3 an ^ / \ 

^ ~ ^ {R^ - Rl)iei~^ - iai~Mn - 1)) J ' 

The estimate 1)2. 47() with this K is valid as long as the denominator in the 
last factor in ()2.53() is > 0, and in order to have a formula for all n we can 
take as a trivial lower bound in other cases or when (|2.47|) is negative. 
As required for ()2.48() . K is monotonously decreasing in n. We now insert 
n = Api^ and obtain 

K{4pf,£) > {I - e) {1 - {l + {const.)Y{£ /af{R^ -Rl)/£^)~^ 

_ (const.)y \ 

{R^ - Rl) {e{a/£f - {const.)Y^{£/a)^) J ^ ' 

with Y = Airpa'^/S as before. Also, the factor 

1 - -L^ = (1 _ {const.)Y'\a/£f) (2.55) 

in (|2.52|) (which is the ratio between n(n — 1) and n^) must not be be 
forgotten. We now make the ansatz 

e ~ y°, a/£ ~ Y'^, {R^ - RI)/£^ ~ (2.56) 

with exponents a, fi and 7 that we choose in an optimal way. The conditions 
to be met are as follows: 

• e{a/£)'^ — (const. )y^(^/a)^ > 0. This holds for all smah enough y, 
provided a + 5/3 < 2 which follows from the conditions below. 

• a > in order that e ^ for y — > 0. 

• 3^ - 1 > in order that Y-'^{a/£f for for y ^ 0. 

• 1 - 3/3 + 7 > in order that Y{£/af{R^-Rl)/£^ for for y ^ 0. 

• 1 — Q — 2/3 — 7>0to control the last factor in (|2.54j) . 
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Taking 

a = 1/17, /3 = 6/17, 7 = 3/17 (2.57) 

all these conditions are satisfied, and 

a = 3/3-l = l-3/3 + 7 = l- a-2/3-7 = 1/17. (2.58) 

It is also clear that 2R/1 ~ Y^l^ = Y^l^"^ ^ up to higher order terms. This 
completes the proof of Theorems 12.31 and 12.41 for the case of potentials 
with finite range. By optimizing the proportionality constants in ()2.56|1 one 
can show that C = 8.9 is possible in Theorem 12.31 ^Sejj. The extension to 
potentials of infinite range but finite scattering length is obtained by approx- 
imation by finite range potentials, controlling the change of the scattering 
length as the cut-off is removed. See Appendix A in |LY2j and Appendix 
B in |LSeYlj for details. We remark that a slower decrease of the potential 
than 1/r^ implies infinite scattering length. ■ 

The exponents ()2.57() mean in particular that 

a < ii « < £ < (2-59) 

whereas Dyson's method required K ^ p^^/^ as already explained. The 
condition p^^l"^ ^ ^ is required in order to have many particles in each box 
and thus n(n — 1) ~ v? . The condition i <C (pa)~^/^ is necessary for a 
spectral gap gap eo(/o) in Temple's inequality. It is also clear that this 
choice of i would lead to a far too big energy and no bound for eo(p) if we 
had chosen Dirichlet instead of Neumann boundary conditions for the cells. 
But with the latter the method works! 



3. The Dilute Bose Gas in 2D 

In contrast to the three-dimensional theory, the two-dimensional Bose gas 
began to receive attention only relatively late. The first derivation of the 
correct asymptotic formula was, to our knowledge, done by Schick [H] for a 
gas of hard discs. He found 

e(p)«i47r^p|ln(/>a2)|-i. (3.1) 

This was accomplished by an infinite summation of 'perturbation series' 
diagrams. Subsequently, a corrected modification of [SI was given in 'HFM' . 
Positive temperature extensions were given in |Poj and in |FHj . All this work 
involved an analysis in momentum space, with the exception of a method 
due to one of us that works directly in configuration space |Llj . Ovchinnikov 
[m derived (jHH)) by using, basically, the method in |Llj . These derivations 
require several unproven assumptions and are not rigorous. 

In two dimensions the scattering length a is defined using the zero energy 
scattering equation (|2.3|) but instead of '(/'('") ~ 1 — a/r we now impose the 
asymptotic condition ^/)(r) ~ ln(r/a). This is explained in the appendix to 
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Note that in two dimensions the ground state energy could not possibly 
be eo(p) ~ 4:TTfipa as in three dimensions because that would be dimen- 
sionally wrong. Since eo{p) should essentially be proportional to p, there is 
apparently no room for an a dependence — which is ridiculous! It turns out 
that this dependence comes about in the ln{pa'^) factor. 

One of the intriguing facts about (|3.1|) is that the energy for N parti- 
cles is not equal to A'^(A'" — l)/2 times the energy for two particles in the 
low density limit — as is the case in three dimensions. The latter quan- 
tity, Eq{2,L), is, asymptotically for large L, equal to 87r//L~^ [ln(L^/a^)] ^ . 
(This is seen in an analogous way as (|2.13j) . The three-dimensional bound- 
ary condition ?/'o(|2?| = R) = 1 — a/R is replaced by V'od^l = R) = ln{R/a) 
and moreover it has to be taken into account that with this normalization 
llV'olP = (volume) (ln(i2/a))^ (to leading order), instead of just the volume in 
the three-dimensional case.) Thus, if the A^(A^— 1)/2 rule were to apply, 1)3. 1(1 
would have to be replaced by the much smaller quantity Airpp [ln(L^/a^)] ^ . 
In other words, L, which tends to oo in the thermodynamic limit, has to 
be replaced by the mean particle separation, p^^/^ in the logarithmic fac- 
tor. Various poetic formulations of this curious fact have been given, but 
the fact remains that the non-linearity is something that does not occur in 
more than two dimensions and its precise nature is hardly obvious, physi- 
cally. This anomaly is the main reason that the two-dimensional case is not 
a trivial extension of the three-dimensional one. 

Eq. (|3.1|) was proved in |LY2j for nonnegative, finite range two-body po- 
tentials by finding upper and lower bounds of the correct form, using similar 
ideas as in the previous section for the three-dimensional case. We discuss 
below the modifications that have to be made in the present two-dimensional 
case. The restriction to finite range can be relaxed as in three dimensions, 
but the restriction to nonnegative v cannot be removed in the current state 
of our methodology. The upper bounds will have relative remainder terms 
0(1 ln{pa^)\^^) while the lower bound will have remainder 0(| ln{pa^)\~^^^). 
It is claimed in ^HFM] that the relative error for a hard core gas is negative 
and 0(ln| ln(/oa^)|| ln(pa^)|~^), which is consistent with our bounds. 

The upper bound is derived in complete analogy with the three dimen- 
sional case. The function /o in the variational ansatz (|2.20() is in two dimen- 
sions also the zero energy scattering solution — but for 2D, of course. The 
result is 

^^/"^ - HbM^'-.pb^ + 0([ln(6/a)]-i)) . (3.2) 

The minimum over b of the leading term is obtained for b = (27r/9)^^/^. 
Inserting this in ()3.2() we thus obtain 

Eo{N, L)/N < j^^^ (1 + 0(1 Hpa'r') . (3.3) 
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To prove the lower bound the essential new step is to modify Dyson's 
lemma for 2D. The 2D version of Lemma 12.51 is: 



Lemma 3.1. Let v{r) > and v{r) = for r > Rq. Let U{r) > he any 
function satisfying 

/•oo 

/ U{r)ln{r/a)rdr < 1 and U{r) = for r < Rq. (3.4) 
Jo 

Let B cM? be star-shaped with respect to (e.g. convex with £ B). Then, 
for all functions ijj in the Soholev space H^{B), 

(/u|VV'(x)p + i?;(|f|)|?/^(x)|2) dx>fi I U{\x\)\ij{x)\^ dx. (3.5) 

Proof. In polar coordinates, r,9, one has jV-i/'p > \^^p/^r\'^. Therefore, it 
suffices to prove that for each angle 6 £ [0,27r), and with ■ilj{r,6) denoted 
simply by /(r), 

R{0) fR{e) 

{fi\df{r)/dr\'^ + lv{r)\f{r)\'^)rdr>fi U{r)\f{r)\'^rdr, (3.6) 



where R{6) denotes the distance of the origin to the boundary of B along 
the ray 6. 

If R{0) < LIq then 1)3. 6|) is trivial because the right side is zero while the 
left side is evidently nonnegative. (Here, u > is used.) 

If R{6) > Rq for some given value of 9, consider the disc ^^{O) = {x £ 
: < |x| < R{6)} centered at the origin in and of radius R{6). Our 
function / defines a spherically symmetric function, x fi\^\) on 1^(9), 
and H3.6|) is equivalent to 

/ f /i|V/(|f |)|2 + ]-v{\x\)\f{\m') dx>J U{\x\)\f{m''dS. 

(3.7) 

Now choose some R G {Rq, R{0)) and note that the left side of ()3.7|) is 
not smaller than the same quantity with 1^(9) replaced by the smaller disc 
= {x £ M? : < |x| < R}. (Again, f > is used.) We now minimize 
this integral over Da, fixing f{R). This minimization problem leads to the 
zero energy scattering equation. Plugging in the solution and integrating by 
parts leads to 

^^V..,..._w..2 , ,,,,,, 



2vr [^^,\df{r)/dr\' + -7;(r)|/(r)|^J rdr > ^-^|/(i?)|^. (3.8) 

The proof is completed by multiplying both sides of 1)3. 8|) by U {R)R ln{R/a) 
and integrating with respect to R from Rq to R{9). ■ 

As in Corollary 12. 6[ Lemma 13.11 can be used to bound the many body 
Hamiltonian Hjsi from below, as follows: 
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Corollary 3.2. For any U as in Lemma \H.l\ and any < e < 1 

Hn > eTn + (1 - e)fiW (3.9) 
with T/v = — /i "^iLi and 



W{xi, . . . ,xn) = y U [ min \xi — Xj\. ] . (3.10) 
For U we choose the fohowing functions, parameterized hy R > Rq: 

UMr) = h''>" f««><'-<^ (3,11) 
I otherwise 

with I'i^R) chosen so that 

/ UR{r)ln{r/a)rdr = 1 (3.12) 
jRo 

for ah R > Rq, i.e., 
fR 

v{R) = / ln(r/a)r dr = \ [R^ {ln{R^/a^) - l) - R^ {ln{Rl/a^) - l) } . 
Jro 

(3.13) 

The nearest neighbor interaction (|3.10|) corresponding to Ur will be denoted 
Wr. 

As in Subsection 12.21 we shall need estimates on the expectation value, 
{Wr)o, of Wr in the ground state of eT/v of ()3.9|1 with Neumann boundary 
conditions. This is just the average value of Wr in a hypercube in M^^. 
Besides the normalization factor v{R), the computation involves the volume 
(area) of the support of Ur, which is 

A{R) = tt{R^ - Rl). (3.14) 

In contrast to the three-dimensional situation the normalization factor 
v{R) is not just a constant {R independent) multiple of A[R); the fac- 
tor ln(r/a) in (|3.4|) accounts for the more complicated expressions in the 
two-dimensional case. Taking into account that Ur is proportional to the 
characteristic function of a disc of radius R with a hole of radius Rq, the 
following inequalities for n particles in a box of side length i are obtained 
by the same geometric reasoning as lead to (|2.36|) . cf. |LYlj : 

{Wr)o > (3-15) 

{Wr)o < (3-16) 
iy[R) L J 

with 

Q = A{R)/e^ (3.17) 
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being the relative volume occupied by the support of the potential Ur. Since 

= v{R)^^Ur we also have 

(W^|)o<-^Wo. (3.18) 



{n - l)Q > fl - (1 - Qt-'A > , I"", (3.19) 



As in [LYT] we estimate [1 - (1 - Q)^""^)] by 

1 + (n - l)Q 
This gives 

n(n — 1) Q , , 

<'^«>» ^ -W- i + («-i)9 - 

n(n — 1) , , 

{Wr)o < \„/ -Q. (3.21) 

From Temple's inequality [Tj we obtain like in (|2.38|) the estimate 

V (»'R>o(ES'-A'(H'fi)ll)/ 

where 

= f (3.23) 
is the energy of the lowest excited state of eT^. This estimate is valid for 
E^^ / > {Wr)o, i.e., it is important that i is not too big. 
Putting ()3.2U|) - H3.22() together we obtain the estimate 

Mn,i)>'^^^K(n) (3.24) 

with 

fl_M)2 / n \ 

= - ^> ■ TTi^ ■ V - (e.(R)IP) - u(„ - 1) J P-^^) 
Note that Q depends on I and ii, and depends on I, R and e besides n. 
We have here dropped the term (Wr)'q in the numerator in (|3.22p . which is 
appropriate for the purpose of a lower bound. 

We note that K is monotonically decreasing in n, so for a given n we may 
replace K{n) by K{p) provided p > n. As explained in the previous section, 
(|2.45|) - (|2.52j) . convexity of n i-^ n{n — 1) together with superadditivity of 
Eq^u, I) in n leads, for p = 4p£^, to an estimate for the energy of N particles 
in the large box when the side length L is an integer multiple of t. 

pA{R)(^ Mw...2^ 



E,{N, L)/N > t-^ Ki^Pn (3.26) 

with p = N/L'^. 

Let us now look at the conditions on the parameters e, R and i that have 
to be met in order to obtain a lower bound with the same leading term as 
the upper bound 
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Prom we have 

We thus see that as long as a < -R < the logarithmic factor in the 

denominator in (|3.27p has the right form for a lower bound. Moreover, for 
Temple's inequality the denominator in the third factor in H3.25() must be 
positive. With n = Apf and iy{R) > (const.) R'^ ln{R'^ /a'^) for R > Rq, this 
condition amounts to 

(const.)eln(i?Va2)/^^ > P^i^- (3.28) 
The relative error terms in (|3.26|) that have to be <C 1 are 

' p^2' ^' ' eRHn{Rya^)- ^ ' 

We now choose 

e~|ln(pa2)|-V5, i ^ p~y^\\n{pa'')\"^'^ , i? ~ p-^/^l ln(pa2)|-Vio 

(3.30) 

Condition ()3.28|) is satisfied since the left side is > (const.) | \n{pa?')\^/^ and 
the right side is ~ | h\{pa'^)\'^f^ . The first three error terms in 1)3. 29() are all 
of the same order, | ln(pa2)|~^/^, the last is ~ | ln(pa^)|~^/^(ln | ln(/>a^)|)~^. 
With these choices, ()3.26|) thus leads to the following: 

Theorem 3.3 (Lower bound). For all N and L large enough such that 
L > (const. )p~^/2| ln(pa2)|Vio and N > (const.)] In(pa2)|^/^ with p = N/L"^, 
the ground state energy with Neumann boundary condition satisfies 

Eo{N, L)/N > ^ij^ (l - 0(1 Hpa'r'/')) . (3.31) 

In combination with the upper bound ()3.3|) this also proves 
Theorem 3.4 (Energy at low density in the thermodynamic limit). 

lim ""^^f , = 1 (3.32) 

pa^^o 4:TTpp\ln[pa^)\-^ 

where eo{p) = lim]\f—,oo Eo{N, p~^^'^N^^'^)/N . This holds irrespective of 
boundary conditions. 



As in the three-dimensional case, Theorem 13.41 is also valid for an infinite 

R 



range potential v provided that f > and for some R we have v{r)r dr < 



oo, which guarantees a finite scattering length. 
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4. BosE-EiNSTEiN Condensation 

Let us now comment on the notion of Bose-Einstein condensation (EEC). 
Given the normahzed ground state wave function ^q{xi, . . . ,xj\f) we can 
form the one-body density matrix which is an operator on L'^{W^) (d = 2 or 
3) given by the kernel 

7(f , x')=N J ^o(x, X)Mx, X)dX , (4.1) 

where we introduced the short hand notation 

N 

X = {x2, . . . , xn) and (iX = J^(ixj. (4.2) 



Then J j{x, x)dx = Tr [7] = A^. BEC in the ground state is the assertion 
that this operator has an eigenvalue of order in the thermodynamic limit. 
Since 7 is a positive kernel and, hopefully, translation invariant in the ther- 
modynamic limit, the eigenfunction belonging to the largest eigenvalue must 
be the constant function L~'^l'^ . Therefore, another way to say that there is 
BEC in the ground state is that 

^ / / ^(x,y)dxdy = 0{X) (4.3) 



as N ^ 00, L ^ 00 with N/L'^ fixed. Unfortunately, this is something that 
is frequently invoked but never proved for many body Hamiltonians with 
genuine interactions — except for one special case: hard core bosons on a 
lattice at half- filling (i.e., N = half the number of lattice sites). The proof 
is in |KLS1. 

The problem remains open after more than 75 years since the first in- 
vestigations on the Bose gas |B| E| • It is also not at all clear that BEC is 
essential for superfluidity, as frequently claimed. Our construction in Sec- 
tionlSlshows that (in 3D) BEC exists on a length scale of order p^^/^y-i/i^ 
which, unfortunately, is not a 'thermodynamic' length like volume^'^'^. The 
same remark applies to the 2D case of Section 3, where BEC is proved over 
a length scale ln(/9a^)|^/^°. 

In a certain dilute limit, however, one can prove (|4.3|) . as has been recently 
shown in |LSej . In this limit the interaction potential v is varied with 
so that the ratio a/L of the scattering length to the box length is of order 
1/A^. In |LSej the case of a Bose gas confined in an external trap potential 
was considered (see Sections El and ^ , but the analysis for a homogeneous 
gas is even simpler and implies the following result. For simplicity, we shall 
treat only the 3D case. 

Theorem 4.1 (BEC in a dilute limit). Assume that, as N ^ 00, p = 
N/L/^ and g = Na/L stay fixed, and impose either periodic or Neumann 
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boundary conditions for H . Then 

lim / / j{x, y)dxdy = 1 . (4.4) 

iv-»oo N L-^ J J 

The reason why the Umit — > oo with Na/L fixed is particularly inter- 
esting will become clear when we study systems confined in a trap potential 
in the next section. Note that the limit we consider is really a limit of a 
dilute gas, since 

as — > oo. Since the ground state energy is of the order ap ^ L^^ is this 
limit, it is also clear why we do not deal with Dirichlet boundary conditions: 
there would be an additional contribution to the energy of the same order, 
and the system would not be homogeneous any more. Dirichlet boundary 
conditions can, however, be treated with the methods of Section |SJ 

At this point we should say what we mean exactly by changing a with N. 
We do this by scaling, i.e., we write 

^(l^l) = 4^i(l^l/a) (4.5) 

for some vi having scattering length 1, and vary a while keeping vi fixed. 
It is easily checked that the v so defined has scattering length a. It is 
important to note that, in the limit considered, a tends to zero (as 
and V becomes a hard potential of short range. This is the opposite of the 
usual mean field limit where the strength of the potential goes to zero while 
its range tends to infinity. 



Proof of Theorem \4-l\ Let g = Na/L and p = N/L^ be fixed as A'' — > oo. 
Since a^p 0, Theorems 12.21 and 12.41 implv 



^lun^ (^^^ J dXdxi U|V^,^'o(^i,X)| 



2 



TV \ 

+i^T;(|fi -fj|)|*o(xi,X)|2 =47r/x5, (4.6) 

j=2 J 

where we again used the short hand notation (|4.2|) . The symmetry of ^'o and 
the boundary conditions have also been used. Even more is true, namely 
that 

lim — dXdxip\Vg^'^o{xi,X)\'^ = Anpgs (4.7) 

N^oo \ p J J 



and 



/Ar\2/3 r. N 
lim^(-) / (iXdxii^i;(|xi -x,|)|^o(^i,^)|' =47r^5(l-5) (4. 
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for some < s < 1. The parameter s is given by s = J | V^'oP/ (4vra), where 
■00 denotes the solution to the scattering equation for v (under the boundary 
condition lim\g\_,^^po{x) = 1; see Eq. (|2.3j) ). This is a simple consequence 
of (|4.6|) by variation with respect to the different components of the energy, 
as was also noted in _CS2. More precisely, the ground state energy is a 
concave function of the mass parameter ^u, so it is legitimate to differentiate 
both sides of H4.6|) with respect to /i. In doing so, it has to be noted that g 
depends on fi through the scattering length. Using (|2.13l) one sees that 

d{fia) 1 



by the Feynman-Hellmann principle, since tpQ minimizes the left side of 

We now concentrate on the term 1)4. 7|) and show that to leading order all 
the energy is located in small balls surrounding each particle. These balls 
can be taken to have radius roughly A^^^/^^ compared to the mean particle 
distance p~^/^ . More precisely, we will show that 

lim A^2/3 ( dX f dfi|V^^^'o(fi,X)|2 = , (4.10) 
where, for fixed X, is given by 



"x 



|x G A min|x -Xfcl > iV""^/^^! . (4.11) 



To see this, we shall show that (|4.6p still holds true if the integrals in 
the first term are restricted to the complement of fi^ , denoted by . The 
proof of this is actually just a detailed examination of the lower bounds to 
the ground state energy derived in Subsection 12.21 What we have to show 
is 

^ j nlVi'^ol^dxidX + j v{\xi-Xj\)\^o\'^dxidX > 4TTfj,apN{l-o{l)) 

(4.12) 



i=l i i<j 

as — > oo, where denotes the set 



n^i = {{xi,X) G A^l mm\xi-Xk\ < N'^/^^} . 

While 1)4.12(1 is not true for all conceivable ^''s satisfying the normalization 
condition ||^'||2 = 1, it is true for ^q. Namely, we claim that 

n\\7i'^\'^dxidX + e / \\7i'^\'^dxidX ] 




XjDl'i'l'^dxidX > 47r^apAf(l - o(l)) (4.13) 



24 



E.H. LIEB, R. SEIRINGER, J. P. SOLOVEJ, AND J. YNGVASON 



for any ^, as long as e > 0{N ^/^'') as ^ oo. Since (|4.7j) implies that 
^'o has total kinetic energy of order 0{apN) = 0{N^/^), (BTT2]) follows from 

It remains to prove (|4.13|) . but this has essentially already been done in 
Subsection l2.2l Namely, the usage of the Dyson Lemma 1231 requires only the 
kinetic energy inside balls of radius R = 0{N^^/^^) around each particle. 
In fact, by Ens. lt^KHll and (ITKTll . R ~ aY'^/^'^ with F ~ ~ and 
a ~ iV-2/3 (ti^g latter because g = Na/L is fixed and L ~ iV^/^). The 
second term in (|4.13|1 . the total kinetic energy multiplied by e, together 
with the 'softened' potential (|2.3fl|) . then gives the desired lower bound to 
the energy (see (|2lSSl)-((23n|), as long as e > ©(A^-^/^^). This proves (|4.13|) 
with an error term of the order 0{N~'^^^'^), and therefore also (|4.12|) with a 
relative error 0{N^'^^^'^). 

Eq. (|4.1()j) means that V£^'o(a^, X) is almost zero (in an sense) outside 
of the small balls $7'^. To conclude BEC we need to show that as a function 

of X, ^oix, X) is essentially constant in $7^. Although has a big volume, 
it can be a weird, and even disconnected, set, so this conclusion is not yet 
possible. However, exploiting the knowledge that the total kinetic energy 
of ^'o (including the balls) is not huge (see (|4.7j) ). we an show the desired 
constancy of ^'o- What we need for this purpose is the following special case 
of the generalized Poincare inequality that will be stated more generally in 
Lemma If). 51 in Section IHl 

Lemma 4.2 (Generalized Poincare inequality, special case). Let A C 

be a cube of side length L, and define the average of a function f G -^^^(A) 

by 

1 



if) A = J3 I fi^) • 



A 



There exists a constant C such that for all measurable sets C A and all 
f G H^{A) the inequality 

I \f{x)-{f)A\'dx<c(L^ [ |V/(f)|2df+|J]f/3 / \Vf{x)\^dx) 
Ja \ Jn J A J 

(4.14) 

holds. Here VL'^ = K\VL, and \ ■ \ denotes the measure of a set. . 

Proof. By scaling, it suffices to consider the case L = 1. Using the usual 
Poincare-Sobolev inequality on A (see |LLoj . Thm. 8.12), we infer that there 
exists a C > such that 

< C^(l|V/lli6/.(^) + ||V/||i./,(^J . (4.15) 
Applying Holder's inequality 

||V/|L«/5(f,) < \\Vfh2^n)M'/^ 
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(and the analogue with 0, replaced by 17^), we see that (|4.14|) holds. ■ 

Applying this result, we are now able to finish the proof of Theorem 14.11 
Denote by (^o)yvx average of ^'o(x, ^) over x £ A. Using Lemma EH 
with = flj^ and f{x) = '^o(x,X) — (\I'o)^j^, we conclude that 

dX I dxho{x,X)-{^o)^xY 



1 dX 






Jn 



(4.16) 



where we used that \^%\ < (47r/3)A^^ 4/i7_ ^Ylie first integral on the right 

^^^^^^ ^^^^^^ 

side of ()4.16|) tends to zero as ^ oo by ()4.1Up . and the second term 
vanishes in this limit because of (|4.7() . Moreover, 



dXdx 



M/o(x,A) 



1 1 



dXdx\mo{x,X)\'' - j^j^ / / 7(f, y)dxdy . (4.17) 



We conclude that 
1 1 



liminf ^ / / ^(x,y)dxdy> lim 

iV— >oo N L'^ J J N-*oo 
and Theorem 14. II is proven. 



dXdx\^o{x,X)\'^ = 1 , (4.18) 



As stated, Theorem 14. II is concerned with a simultaneous thermodynamic 
and a — > limit, where, as A^ ^ oo, the box length L is proportional to N^^^ 
and a ~ A^~^/^. By scaling, the above result is equivalent to considering 
a Bose gas in a fixed box of side length L = 1, and keeping Na fixed as 
A^ — > oo, i.e., a ~ 1/A^. The ground state energy of the system is then, 
asymptotically, A^ x AirNa, and Theorem 14.11 implies that the one-particle 
reduced density matrix 7 of the ground state converges, after division by 
A^, to the projection onto the constant function. An analogous result holds 
true for inhomogeneous systems. This was recently shown in |LSej and will 
be presented in Sectional 



5. Gross-Pitaevskh Equation for Trapped Bosons 

In the recent experiments on Bose condensation (see, e.g., |KDj ) . the par- 
ticles are confined at very low temperatures in a 'trap' where the particle 
density is inhomogeneous^ contrary to the case of a large 'box', where the 
density is essentially uniform. We model the trap by a slowly varying con- 
fining potential V , with V{x) — > 00 as \x\ — > 00. The Hamiltonian becomes 
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N 

H = ^{-fiAi + V{xi)}+ v{\xi-Xj\). (5.1) 

i=l l<i<j<N 

Shifting the energy scale if necessary we can assume that V is nonnegative. 
The ground state energy, f\w, of —//A + V{x) is a natural energy unit and 
the corresponding length unit, ^Jh/ {muj) = \j2iil(fuxi) = Lose, is a measure 
of the extension of the trap. 

In the sequel we shall be considering a limit where a/Losc tends to zero 
while N ^ oo. Experimentally a/Losc can be changed in two ways: One 
can either vary Lose ^* The first alternative is usually simpler in prac- 
tice but very recently a direct tuning of the scattering length itself has 
also been shown to be feasible [CCRCWj . Mathematically, both alterna- 
tives are equivalent, of course. The first corresponds to writing V{x) = 
-^osc^i(^/-^osc) and keeping Vi and v fixed. The second corresponds to writ- 
ing the interaction potential as v(|x|) = a^'^vi{\x\/a) like in (|4.5|) . where 
vi has unit scattering length, and keeping V and vi fixed. This is equiva- 
lent to the first, since for given Vi and vi the ground state energy of (|5.1|) . 
measured in units of huj, depends only on N and a/Losc- In the dilute limit 
when a is much smaller than the mean particle distance, the energy becomes 
independent of vi. 

We choose Lose as a length unit. The energy unit is hio = 2ij,L^^c = 
Moreover, we find it convenient to regard V and vi as fixed. This justifies 
the notion Eq (N, a) for the quantum mechanical ground state energy. 

The idea is now to use the information about the thermodynamic limiting 
energy of the dilute Bose gas in a box to find the ground state energy of 
(|5.1j) in an appropriate limit. This has been done in |LSeYH |LSeY2j and in 
this section we give an account of this work. As we saw in Sections |2l and IHl 
there is a difference in the p dependence between two and three dimensions, 
so we can expect a related difference now. We discuss 3D first. 

5.1. Three Dimensions. Associated with the quantum mechanical ground 
state energy problem is the Gross-Pitaevskii (GP) energy functional |G1UG21 

El 

£G^[cj)]= [ {fi\V<P\^ + V\<P\^ + i7rfia\<P\^)dx (5.2) 
with the subsidiary condition 

/ \^\' = N. (5.3) 

As before, a > is the scattering length of v. The corresponding energy is 
i?GP(Ar,a)= inf fGP[0]=f«P[0GP], (5.4) 

with a unique, positive c/)^^. The existence of the minimizer (f>'^^ is proved 
by standard techniques and it can be shown to be continuously differentiable. 
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see |LSeYlj . Sect. 2 and Appendix A. The minimizer depends on N and a, 
of course, and when this is important we denote it by i^^^. 

The variational equation satisfied by the minimizer is the GP equation 

-^A0GP(f) + V{x)(t>^^{x) + 8Trfiacl)^P{xf = fi^^ c/)^^ (x) , (5.5) 

where is the chemical potential, given by 

11^^ = dE^^{N,a)/dN = E^^{N,a)/N + {4TTf,a/N) [ {x)\Ux. (5.6) 



The GP theory has the following scaling property: 

E^^{N, a) = NE^^{1, Na), (5.7) 

and 

^Z(^) = N'/'4>?S^,{x). (5.8) 
Hence we see that the relevant parameter in GP theory is the combination 
Na. 

We now turn to the relation of E'^^ and c/)'^^ to the quantum mechanical 
ground state. If u = 0, then the ground state of 1)5. 1(1 is 

^'o(xi,... ,xn) = l\f=i<Po{xi) 

with (j)Q the normalized ground state of — /iA + V{x). In this case clearly 
= \/iV (f)Q, and then E'^^ = Ntiw = Eq. In the other extreme, if 
V{x) = for X inside a large box of volume L'^ and V{x) = oo otherwise, 
then PS ^N/L^ and we get E'^^{N,a) = An jiaN'^ / , which is the 
previously considered energy Eo for the homogeneous gas in the low density 
regime. (In this case, the gradient term in S'^^ plays no role.) 
In general, we expect that for dilute gases in a suitable limit 

Eq^E"^^ and pQ^(f)p« |0GP(f)|' = pGP(f), (5.9) 

where the quantum mechanical particle density in the ground state is defined 
by 

p^^{x)=N J \^o{x,X2,...,XN)\'^dX2---dXN- (5.10) 

Dilute means here that 

pa^ < 1, (5.11) 

where 

P=^l \p''''m'dx (5.12) 

is the mean density. 

The limit in which (|5.9j) can be expected to be true should be chosen so 
that all three terms in S'^^ make a contribution. The scaling relations (|5.7|) 
and (|5.8j) indicate that fixing Na as ^ oo is the right thing to do (and this 
is quite relevant since experimentally A^ can be quite large, 10^ and more, 
and Na can range from about 1 to 10^ (DGPS ). Fixing Na (which we refer 
to as the GP case) also means that we really are dealing with a dilute limit. 
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because the mean density p is then of the order N (since pN,a = ^Pi,Na) 
and hence 

a^p ~ (5.13) 
The precise statement of (|5.9|) is: 

Theorem 5.1 (GP limit of the QM ground state energy and den- 
sity). If N ^ oo with Na fixed, then 

hm "\ ' \ = 1, (5.14) 

and 

hm ^pZi^) = \<P?S^a(x)\' (5.15) 

iv— >oo iV 

in the weak -sense. 

To describe situations where Na is very large, it is appropriate to consider 
a Hmit where, as — > oo, a ^ iV^^, i.e. Na oo, but still pa'^ — > 0. In this 
case, the gradient term in the GP functional becomes negligible compared 
to the other terms and the so-called Thomas-Fermi (TF) functional 

£^^[p]= [ {Vp + 4TTpap'^)dx (5.16) 

arises. (Note that this functional has nothing to do with the fermionic 
theory invented by Thomas and Fermi in 1927, except for a certain formal 
analogy.) It is defined for nonnegative functions p on M'^. Its ground state 
energy E'^^ and density p^^ are defined analogously to the GP case. (The 
TF functional is especially relevant for the two-dimensional Bose gas. There 
a has to decrease exponentially with in the GP limit, so the TF limit is 
more adequate; see Subsection 15 . 21 below) . 

Our second main result of this section is that minimization of ()5.16|) re- 
produces correctly the ground state energy and density of the many-body 
Hamiltonian in the limit when N — > oo, a^p — > 0, but Na oo (which we 
refer to as the TF case), provided the external potential is reasonably well 
behaved. We will assume that V is asymptotically equal to some function 
W that is homogeneous of some order s > 0, i.e., W{Xx) = X'^W{x) for all 
A > 0, and locally Holder continuous (see |LSeY2j for a precise definition). 
This condition can be relaxed, but it seems adequate for most practical 
applications and simplifies things considerably. 

Theorem 5.2 (TF limit of the QM ground state energy and den- 
sity). Assume that V satisfies the conditions stated above. If g = Na oo 
as N ^ oo, but still a^p — > 0, then 

:^'^oo E^^{N,a) =^^ ^'-'^^ 



and 



hm '-^pTaia'^^'"-'^^) = pJJiS) (5.18) 

iV— >oo iV ' 
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in the weak L -sense, where pli is the minimizer of the TF functional under 
the condition J p = 1, a = 1, and with V replaced by W . 

In the following, we will present the essentials of the proofs Theorems 15.11 
and 15.21 We will derive appropriate upper and lower bounds on the ground 
state energy Eq. 

The proof of the lower bound in Theorem 15 . II presented here is a modified 
version of (and partly simpler than) the original proof in [LSeYlj . 

The convergence of the densities follows from the convergence of the en- 
ergies in the usual way by variation with respect to the external potential. 
For simplicity, we set ^ = 1 in the following. 

Proof of Theorem,s \5.1\ and \5.Si Part 1: Upper hound to the QM energy. To 
derive an upper bound on £^0 we use a generalization of a trial wave function 
of Dyson jPlj . who used this function to give an upper bound on the ground 
state energy of the homogeneous hard core Bose gas (c.f. Section f2.1jl . It is 
of the form 

N 

^{xi, . . . ,xn) =J{<f^ {x^)F{xi, . . . ,xn), (5.19) 

i=l 

where F is constructed in the following way: 

N 

F{xi, . . . ,xn) = J{f{ti{xi, . . . ,Xi)), (5.20) 

i=l 

where ti = min{|xj — Xj\,l < j < i — 1} \s the distance of xi to its nearest 
neighbor among the points xi, . . . , and / is a function of r > 0. As in 
1)2. 2U() we choose it to be 

/(r) = I /o(0//o(6) for r <b 

•I "-^ ' Y 1 for r > 6, ^ ' 

where /o is the solution of the zero energy scattering equation (|2.l'{j) and b 
is some cut-off parameter of order b ~ p^^^^. The function (|5.19j) is not 
totally symmetric, but for an upper bound it is nevertheless an acceptable 
test wave function since the bosonic ground state energy is equal to the 
absolute ground state energy. 

The result of a somewhat lengthy computation (see |LSeYlj for details) 
is the upper bound 

Eq{N, a) < E'^^iN, a) (l + 0{ap^/^)^ . (5.22) 

Part 2: Lower bound to the QM energy, GP case. To obtain a lower bound 
for the QM ground state energy the strategy is to divide space into boxes and 
use the estimate on the homogeneous gas, given in Theorem 1231 in each box 
with Neumann boundary conditions. One then minimizes over all possible 
divisions of the particles among the different boxes. This gives a lower bound 
to the energy because discontinuous wave functions for the quadratic form 
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defined by the Hamiltonian are now allowed. We can neglect interactions 
among particles in different boxes because f > 0. Finally, one lets the 
box size tend to zero. However, it is not possible to simply approximate 
V hy a constant potential in each box. To see this consider the case of 
noninteracting particles, i.e., v = and hence a = 0. Here Eq = NtkJ, but 
a 'naive' box method gives only mmgV{x) as lower bound, since it clearly 
pays to put all the particles with a constant wave function in the box with 
the lowest value of V. 

For this reason we start by separating out the GP wave function in each 
variable and write a general wave function ^ as 

N 

^{xi, . . . ,xn) = Yl(^^^{xi)F{xi, . . . ,xn). (5.23) 

i=l 

Here cp'^^ = c/)^^ is normalized so that / \(p'^^\'^ = N. Eq. defines F 

for a given ^' because (p^^ is everywhere strictly positive, being the ground 
state of the operator —A + V + 87ra\(l)^^\'^. We now compute the expectation 
value of H in the state ^. Using partial integration and the variational 
equation (|5.5|) for cp'^^ , we see that 



with 

N 

Q{F) = 



^''''(^'«) = 4™ j |p«P|2 + Q(F), (5.24) 



(5.25) 

We recall that p'^^{x) = \4>Wai^)\'^- ^'^^ computing the ground state energy 
of H we have to minimize the normalized quadratic form Q. Compared to 
the expression for the energy involving ^ itself we have thus obtained the 
replacements 

N N 

V{x) -87rap^P(x) and JJtfxi ^ Ylp'^'^ {xi)dxi . (5.26) 

i=l i=l 

We now use the box method on this problem. More precisely, labeling the 
boxes by an index a, we have 



infQ(F) > inf Vinf g„(F„), (5.27) 

F {na } — ^ FoL 

a 

where Qa is defined by the same formula as Q but with the integrations 
limited to the box a, is a wave function with particle number tIq,, and 
the infimum is taken over all distributions of the particles with ^ = A^. 

We now fix some M > 0, that will eventually tend to oo, and restrict 
ourselves to boxes inside a cube A^f of side length M. Since u > the 
contribution to (|5.27jl of boxes outside this cube is easily estimated from 
below by — STrA^asup^^y^M P'^^ i^)^ which, divided by A^, is arbitrarily small 



(5.30) 



THE GROUND STATE OF THE BOSE GAS 31 

for M large, since Na is fixed and (jP^ /N^/'^ = 4>i^^ decreases faster than 
exponentially at infinity ( |LSeYlj . Lemma A. 5). 

For the boxes inside the cube Km we want to use Lemma lTKl and therefore 
we must approximate p^^ by constants in each box. Let Pa,max and Pa^min-, 
respectively, denote the maximal and minimal values of p^^ in box a. Define 

ria 

. . . ,x„J = . . . ,x„J JJ 4>^^{xk), (5.28) 

k=l 

and 

Tie, 

^'«(fi,...,f„J =F„(xi,...,x„J J]0GP(xfc). (5.29) 

fc=i 

We have, for all 1 < i < tIq, 

/ n P'^'^ixk) (|V.F„|2 + 1 ^(|f, - f,|)|F«|2) 

k=l j^i 

We now use Lemma 12.51 to get, for all < e < 1, 

m^>PaMnj (£|V^^'«P + a(l-e)C/(ti)|*L'^P) (5.31) 

where U is the distance to the nearest neighbor of Xi, c.f., ()2.29|) . and U the 
potential ^E^i . 

Since ^'^ = 4>'^^ {xi)'^a^ we can estimate 

IV.^-^P < 2p„,^ax|V,^«|2 + 2|*»|2iVCM (5.32) 

with 

Cm = ^ sup |V</.GP(x)|2 = sup \V4>?!n,{x)\\ (5.33) 

Since Na is fixed. Cm is independent of A^. Inserting (|5.32p into ()5.3ip . 
summing over i and using /9^^(xj) < Pa.max in the last term of (|5.25j) (in 
the box a), we get 

Qa{Fa) > ^"'""" -gf (n-a, L) - 87rapa,maxna - sCMna, (5.34) 
Po,max 

where L is the side length of the box and E^{na,L) is the ground state 
energy of 

^i-^eAi + {l-e)aUiU)) (5.35) 
j=i 

in the box (c.f. (|2.35|) ). We want to minimize (|5.34j) with respect to and 
drop the subsidiary condition = A in (|5.27l) . This can only lower the 

minimum. For the time being we also ignore the last term in (|5.34)1 . (The 
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total contribution of this term for all boxes is bounded by eCmN and will 
be shown to be negligible compared to the other terms.) 

Since the lower bound for the energy of Theorem 12.41 was obtained pre- 
cisely from a lower bound to the operator (|5.35|) . we can use the statement 
and proof of Theorem 12.41 From this we see that 



£;fK,L) > (l-.)^(l-CyiA^) (5.36) 
with Ya = a?na/L^, provided Ya is small enough, e > Ya'^'^ and 

— 1/17 

(const. )yj3- . The condition on e is certainly fulfilled if we choose e = 
yi/i'^ with Y = a^N/L^ . We now want to show that the no minimizing the 
right side of (|5.34|) is large enough for (|5.3(ij) to apply. 

If the minimum of the right side of (|5.,'-{4)1 (without the last term) is taken 
for some n„, we have 



P«,min / j-^U 



E^ifia + l,L)- E^{na,L)) > 87rapa,max. (5.37) 



PQ,max 

On the other hand, we claim that 
Lemma 5.3. For any n 

E^{n + 1,L)- E^{n,L) < Svra-^. (5.38) 

Proof. Denote the operator 1)5. 35(1 by Hn, with Ua = n, and let be its 
ground state. Let t'^ be the distance to the nearest neighbor of among 
the n + 1 points xi, . . . , Xn+i (without Xi) and ti the corresponding distance 
excluding Xn+i- Obviously, for 1 < z < n, 



and 



Therefore 



U{t'i)<U{U) + U{\xi.-Xn+i\) (5.39) 

n 

U{t',,+,)<Y,U{\xi-Xn+i\). (5.40) 



i=l 



Hn+1 <Hn- + 2a^ [/(|xi - Xn+l\). (5.41) 

i=l 

Using ^'n/L'^/^ as trial function for Hn+i we arrive at ()5.38|) . ■ 
Eq. (|08|) together with ((07|) shows that TJq, is at least ~ Pa,max 

L3. We 

shall choose L ~ A^~^/^'^, so the conditions needed for ()5.36|) are fulfilled for 
A'' large enough, since Pa,ia&x ~ N and hence fia ~ N'^/^^ and Ya ~ A^~^. 
In order to obtain a lower bound on Qa we therefore have to minimize 

4vra ( { 1 - cyVir) _ Sn^p^,^,,") . (5.42) 



Pa, max 
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We can drop the requirement that has to be an integer. The minimum 
of (15.421) is obtained for 



^ Pa, max ^^-^ 

Pa, mill (1 CY^I^'^^ 

By Eq. 1)5. 24() this gives the following lower bound, including now the last 
term in (|5.34|) as well as the contributions from the boxes outside Km, 

Eo{N,a) - E^^{N,a) > 

aCAAf \Pa,minV^ 

- y^/i^iVCAf - 47raiV sup p^^(f). 

Now yO*-'^ is differentiable and strictly positive. Since all the boxes are in the 
fixed cube Am there are constants C < oo, C" > 0, such that 

/5a, max ~ Pa,min ^ NC L, Pa,iain ^ NC . (5.45) 

Since L ~ A^^^/^o and Y ~ N^^"^^^^ we therefore have, for large A^, 

i=a3c^^^-(--)iV-V» (5.46) 

Also, 

4™ ^ /.^ „^i,L3 < ina j jp^^l^ < E''''{N,a). (5.47) 

acAM 

Hence, noting that E'^^{N, a) = NE'^^{1, No) ~ N since Na is fixed, 

i^Gpf^"\ > 1 - (const.)(l + CM)iV-^/^° - (const.) sup \ct>%f, (5.48) 

where the constants depend on Na. We can now take N ^ oo and then 
M ^ oo. 

Pari 5.' Lower bound to the QM energy, TF case. In the above proof of 
the lower bound in the GP case we did not attempt to keep track of the 
dependence of the constants on Na. In the TF case Na ^ oo, so one would 
need to take a closer look at this dependence if one wanted to carry the 
proof directly over to this case. But we don't have to do so, because there is 
a simpler direct proof. Using the explicit form of the TF minimizer, namely 

pZi^) = ^i^^^^-yi^)]+^ (5.49) 

where = max{t, 0} and fx^^ is chosen so that the normalization condi- 
tion f pj^^ = N holds, we can use 

V{x) > ^i'^^ - 8TTap^^{x) (5.50) 
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to get a replacement as in (|5.26j) . but without changing the measure. More- 
over, p^^ has compact support, so, applying again the box method described 
above, the boxes far out do not contribute to the energy. However, fi^^ 
(which depends only on the combination Na) tends to infinity as Na — > oo. 
We need to control the asymptotic behavior of f^^^, and this leads to the 
restrictions on V described in the paragraph preceding Theorem 15.21 For 
simplicity, we shall here only consider the case when V itself is homogeneous, 
i.e., V{Xx) = X^V{x) for all A > with some s > 0. 
In the same way as in (|5.()|) we have, with g = Na, 



fi''ig)=dE''iN,a)/dN = E''il,g)+47rg J Ipllix^dx. (5.51) 
The TF energy, chemical potential and minimizer satisfy the scaling relations 

E^^il,g)=g-^/(^-^'^E^^il,l), (5.52) 

t^^^ig) = //(^+='V^"(i), (5.53) 

and 

^3/(s+3)^TF(^l/(s+3)^^) = Pr,{S). (5.54) 

We also introduce the scaled interaction potential, v, by 

v{x) = 52/(^+3)^;(5i/(«+3)^) (5_55) 

with scattering length 

2 = 5-i/{^+3)a. (5.56) 
Using (|5.5()|) . (|5.51|) and the scaling relations we obtain 

Eo{N, a) > i?TF(Ar, a) + 47riV5^/(^+3) J \pJJ\^ + g-^^'+^^Q (5.57) 

with 

Q = inf_ [ d^i^l" +hYl ^(^^ - - 8nNapJJ{xi)\^\^) . 

(5.58) 

We can now proceed exactly as in Part 2 to arrive at the the analogy of Eq. 



.44|) . which in the present case becomes 
Eo{N,a) - E'^^{N,a) > 

Here Pa,ma.x is the maximum of pj^ in the box a, and Y = a^N/L^. This 
holds as long as L does not decrease too fast with N. In particular, if L is 
simply fixed, this holds for all large enough N. Note that 

p = Npi,g^Ng-^/^'+^^pi,u (5.60) 
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SO that a^N ~ a^p goes to zero as — > oo by assumption. Hence, if we first 
let — > oo (which imphes y — > 0) and then take L to zero, we of arrive at 
the desired result 

hminf >1 (5.61) 

in the limit — > oo, a^p — > 0. Here we used the fact that (because V, and 
hence p^^, is continuous by assumption) the Riemann sum X^o Pa.max-^^ 
converges to /|Pi"TP as L ^ 0. Together with the upper bound 1)5. 22() 
and the fact that E^^ {N , a) / E^^ {N , a) = E^^{l,Na)/E^^{l,Na) 1 as 
Na oo, which holds under our regularity assumption on V (c.f. Lemma 
2.3 in |LSeY2p . this proves and (IKTTj) . 

Part 4-' Convergence of the densities. The convergence of the energies im- 
plies the convergence of the densities in the usual way by variation of the 
external potential. We show here the TF case, the GP case goes analogously. 
Set again g = Na. Making the replacement 

V{x) V{x) + 55^/(^+3)z((7-i/(^+3)x) (5.62) 

for some positive Z £ and redoing the upper and lower bounds we see 
that (|5.17|) holds with W replaced hy W + SZ. Differentiating with respect 
to 5 at (5 = yields 

hm '^^pZia'^''"-'^^) = PiJ(^) (5.63) 

in the sense of distributions. Since the functions all have L^-norm 1, we can 
conclude that there is even weak -convergence. ■ 

5.2. Two Dimensions. In contrast to the three-dimensional case the en- 
ergy per particle for a dilute gas in two dimensions is nonlinear in p. In view 
of Schick's formula (|3.H) for the energy of the homogeneous gas it would ap- 
pear natural to take the interaction into account in two dimensional GP 
theory by a term 

4tt [ |ln(|(/.(x)|2a2)|-i|0(f)|4df, (5.64) 
and such a term has, indeed, been suggested in |Shj and KJNSCj . However, 



since the nonlinearity appears only in a logarithm, this term is unnecessarily 
complicated as far as leading order computations are concerned. For dilute 
gases it turns out to be sufficient, to leading order, to use an interaction 
term of the same form as in the three-dimensional case, i.e, define the GP 
functional as (for simplicity we put p = 1 in this section) 

£^^[(t>]= ! (|V(/>|2 + y|(/)|2 + 47ra|(/>|^)dx, (5.65) 

where instead of a the coupling constant is now 

a = \\a{pNa^)\~^ (5.66) 
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with the mean density for the GP functional at coupling constant 1 and 
particle number A^. This is defined analogously to (|5.12j) as 

~PN = l^j \&dS (5-67) 

where is the minimizer of (|5.65|) with a = \ and subsidiary condition 
/|(/>P = A^. Note that a in (|5.66|) depends on through the mean density. 

Let us denote the GP energy for a given N and coupling constant a by 
E'^^{N, a) and the corresponding minimizer by i;^^^. As in three dimensions 
the scaling relations 

E^^{N,a) = NE^^{l,Na) (5.68) 

and 

hold, and the relevant parameter is 

g = Na. (5.70) 

In three dimensions, where a = a, it is natural to consider the limit — > 
oo with g = Na= const. The analogue of Theorem l5.1l in two dimensions is 

Theorem 5.4 (Two-dimensional GP limit theorem). If, for N oo, 

o-^pN with g = N/\ h\{a?pi~j)\ fixed, then 

Eo{N,a) _ 
N-^oo ^GP(^r,l/|ln(a2p^)|) ~ ^ ^ ' 



and 



in the weak -sense. 



This result, however, is of rather limited use in practice. The reason is 
that in two dimensions the scattering length has to decrease exponentially 
with N if g is fixed. The parameter g is typically very large in two dimensions 
so it is more appropriate to consider the limit N ^ oo and g ^ oo (but still 
pNa^ 0). 

For potentials V that are homogeneous functions of x, i.e., 

ViXx) = X'Vix) (5.73) 

for some s > 0, this limit can be described by the a 'Thomas-Fermi' energy 
functional like (|5.16|) with coupling constant unity: 

£^^[p]= [ {V{x)p{x) + 4TTp{xf)dx. (5.74) 

This is just the GP functional without the gradient term and a = 1. Here 
p is a nonnegative function on and the normalization condition is 

p{x)dx = 1. (5.75) 
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The minimizer of (|5.74j) can be given explicitly. It is 



■^^^^-^-{8nrV'-V{x)U (5.76) 



where the chemical potential fi^^ is determined by the normalization con- 
dition (|5.75|) and = t for t > and zero otherwise. We denote the 
corresponding energy by E'^^il,!). By scaling one obtains 

hm ii;«P(l,5)/<?^/(^+') = ^^^(1, 1), (5.77) 

lim 9'/^'^'^P?^,{9'/^-'^'^x) = pJl{x), (5.78) 

with the latter limit in the strong L? sense. 

Our main result about two-dimensional Bose gases in external potentials 
satisfying (|5.7Hj) is that analogous limits also hold for the many-particle 
quantum mechanical ground state at low densities: 

Theorem 5.5 (Two-dimensional TF limit theorem). In two dimen- 
sions, if a^pM — > 0, but g = N/\ ln(/>7vfl^)| — > cxd as N ^ oo then 

1- ^^ = ^""(M) (5.79) 

N—>00 g^/^^^ 

and, in the weak sense, 

lim '^^pZia'^^'"-'^-) = PJJi-)- (5-80) 

Remarks: 1. As in Theorem 15.21 it is sufficient that V is asymptotically 
equal to some homogeneous potential, W. In this case, S"'"^(l,l) and pj^ 
in Theorem 15.51 should be replaced by the corresponding quantities for W. 
2. From Eq. it follows that 

PN ~ iV^/(^+2) (5.81) 

for large A^. Hence the low density criterion a'^p <^ 1, means that a/ Lose ^ 

We shall now comment briefly on the proofs of Theorems 15.41 and 15.51 
mainly pointing out the differences from the 3D case considered previously. 

The upper bounds for the energy are obtained exactly in a same way as in 
three dimensions. For the lower bound in Theorem 15.41 the point to notice 
is that the expression (|5.42|) . that has to be minimized over n^, is in 2D 
replaced by 

1 Pa,ra\n ^oc ^ f ^ \ '^^aPo 



L2 I In(a2n<,/L2)| | ln{a'^N/ L^)]^^ J \ Ha^pN) 

(5.82) 

since Eq. (|5.36|) has to be replaced by the analogous inequality for 2D (c.f. 
(|3.31|) ). To minimize (|5.82|) we use the following lemma: 
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Lemma 5.6. For < x, 6 < 1 and k > 1 we have 

6 . , 6^ 1 



, 2——xk > -—— 1 + , , , , ■ (5.83) 

\lnx\ \lnb\ - \lnb\ \ (2|ln6|)2 J ^ ' 

Proof. Replacing x hy x/k and using the monotonicity of In we see that it 
suffices to consider k = 1. Since Inx > —^x^'^ for all d > we have 

_ 2^ > \}^edx^+'' - ^ > c(d)(6'^ed| ln6|)-i/(i+'^) (5.84) 
0^ I lnx| 0^ b 

with 

cid) = 2(2+'^)/(i+^) ( I i ^ > -1 - -d^ 

^ ' \i2 + d)(^+<i)/i^+d) (2 + d)i/(i+d)y - 4 

(5.85) 

Choosing d = 1/| ln5| gives the desired result. ■ 
Applying this lemma with x = a'^ria/L'^, b = a'^pa,max and 

Pa,miix A ^ \ I ln(a^Po max)| 



Pa.min V | ln(a2Ar/L2) 1 1/5 ^ |ln(a2pjv)| 
we get the bound 

n2 T'^ 



(EH > -47r ,;™_ ,| 1 + ^^-^ fc. (5.87) 

|ln(a2p^)| V 4|ln(a2p„,j„ax)IV 

In the limit considered, k and the factor in parenthesis both tend to 1 and 
the Riemann sum over the boxes a converges to the integral as L — > 0. 

The TF case, Thm. l?3l is treated in the same way as in three dimensions, 
with modifications analogous to those just discussed when passing from 3D 
to 2D in GP theory. 



6. BEG FOR Dilute Trapped Gases 

It is gratifying to see the experimental realization, in traps, of the long- 
predicted Bose-Einstein condensation (BEG) of gases. From the theoretical 
point of view, however, a rigorous demonstration of this phenomenon - 
starting from the many-body Hamiltonian of interacting particles - has not 
yet been achieved. Following [LSij, we will provide in this section such a 
rigorous justification for the ground state of 2D or 3D bosons in a trap 
with repulsive pair potentials, and in the well-defined limit in which the 
Gross-Pitaevskii (GP) formula is applicable. It is the first proof of BEG for 
interacting particles in a continuum (as distinct from lattice) model and in 
a physically realistic situation. The Gross-Pitaevskii limit under discussion 
here is, of course, a physically simpler limit than the usual thermodynamic 
limit in which the average density is held fixed as the particle number goes 
to infinity. In the GP limit one also lets the range of the potential go to zero 
as goes to infinity, but in such a way that the overall effect is non-trivial. 
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That is, the combined effect of the infinite particle hmit and the zero range 
limit is such as to leave a measurable residue — the GP function. 

It was shown in the previous section (see also Theorem l6.3l below) that, for 
each fixed Na, the minimization of the GP functional correctly reproduces 
the large asymptotics of the ground state energy and density oi H - 
but no assertion about BEG in this limit was made. We will now extend 
this result by showing that in the Gross-Pitaevskii limit there is indeed 
100% Bose condensation in the ground state. This is a generalization of 
the homogeneous case considered in Theorem 14.11 In the following, we 
concentrate on the 3D case, but analogous considerations apply also to the 
2D case. 

For use later, we define the projector 

pGP^|^GP^^^GP| _ (g_i) 

Here (and everywhere else in this section) we denote ^p^^ = (f^i^a ^"-"^ sim- 
plicity, where (/'j^^a minimizer of the GP functional ()5.2() with pa- 
rameter Na and normalization condition J = 1 (compare with (^IHI))- 
Moreover, we set /i = 1. 

In the following, \I'o denotes the (nonnegative and normalized) ground 
state of the Hamiltonian (|5.1() . BEG refers to the reduced one-particle den- 
sity matrix 7(x, x') of ^'o, defined in (|4.1|) . 

Gomplete (or 100%) BEG is defined to be the property that ;^7(x, x') not 
only has an eigenvalue of order one, as in the general case of an incomplete 
BEG, but in the limit it has only one nonzero eigenvalue (namely 1). Thus, 
-^7(x, x') becomes a simple product Lp{x)*Lp{x') as — > oo, in which case 
(p is called the condensate wave function. In the GP limit, i.e., — > oo 
with Na fixed, we can show that this is the case, and the condensate wave 
function is, in fact, the GP minimizer ipP^ . 

Theorem 6.1 (Bose-Einstein condensation). For each fixed N a 
lim ^Mx,x') = 99^P(x)v;^P(f') . 

in trace norm, i.e., Tr |-^7 — 0. 

We remark that Theorem implies that there is also 100% condensation 
for all n-particle reduced density matrices 

(n) / -' . \ 

7 \X\, . . . , Xn\ X^, . . . , X^j 

^ / ^o(^l,---,^Af)*o(5l,---,^n;^n+l,---^Af)d^n+l---rf^Ar 

(6.2) 

of ^'O) i-e., they converge, after division by the normalization factor, to the 
one-dimensional projector onto the n-fold tensor product of (/?^^. In other 
words, for n fixed particles the probability of finding them all in the same 
state ip'^^ tends to 1 in the limit considered. To see this, let a* ,a denote the 
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boson creation and annihilation operators for the state and observe 

that 

1 > hm Af-"(^o|(«*)"a"|^o) = lim iV-"(^'o|(a*a)"|^'o) , (6.3) 

N^oo N—*oo 

since the terms coming from the commutators [a, a*] = 1 are of lower order 
as — > oo and vanish in the limit. From convexity it follows that 

iV-"(^'o|(a*a)"|^'o) > iV""(^o|a*a|^o>" (6.4) 

which converges to 1 as — > oo, proving our claim. 

Another corollary, important for the interpretation of experiments, con- 
cerns the momentum distribution of the ground state. 

Corollary 6.2 (Convergence of momentum distribution). Let 

p(k) = J J ^{x, x) exp[ik ■ (x — x)]dxdx 
denote the one-particle momentum density 0/^*0 • Then, for fixed Na, 

hm Ip(fc) = |^«P(fc)P 

strongly in L^(M^). Here, denotes the Fourier transform of cp'^^ . 

Proof. If J-' denotes the (unitary) operator 'Fourier transform' and if h is an 
arbitrary L°°-function, then 



1 



-- \TT[J^-\j/N-P^^)J^h]\ 
< ||/i|UTV|7/Ar-pGP|^ 

from which we conclude that 

||p/iV- < Tr|7/Ar-pGP| _ 



Before proving Theorem 16.11 let us state some prior results on which we 
shall build. Then we shall formulate two lemmas, which will allow us to 
prove Theorem 16. II 

The following theorem is an extension of Theorem 15. II 

Theorem 6.3 (Asymptotics of the energy components). If ipo denotes 
the solution to the zero-energy scattering equation for v (under the boundary 
condition lim|^|_>oo ■0o(^) = and s = J |VV'o|^/(47ra), then < s < 1 and 
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lim / \Vg^'^o{xi,X)\^dxidX 

■f^oo J 

= J \V^^^{x)\^dx + ATrNas j \ip^^ {x)\'^dx, (6.5a) 

lim I V{xi)\^o{xi,X)\^dxidX = f V {x)\^^^ {x)\^ dx , (6.5b) 

lim / v{\xi - Xj\)\'iio{xi,X)\'^dxidX 



i=2 



= (1 - s)iTTNa J \ip^^{x)\^dx. (6.5c) 
Here we introduced again the short hand notation 

N 

X = {x2, ■ ■ ■ , xj\f) and (iX = J^dxj. (6-6) 

i=2 

Theorem 16.31 is a simple consequence of Theorem 15.11 by variation with 
respect to the different components, in the same way as was explained in 
the proof of Theorem lO c.f. Eqs. (|I?7|) - (|0)) . This was also noted in |CS2j . 

As already stated, Theorem 16.11 is a generalization of Theorem 14.11 the 
latter corresponding to the case that F is a box potential. It should be noted, 
however, that we use different scaling conventions in these two theorems: In 
Theorem 14. II the box size grows as N^^^ to keep the density fixed, while in 
Theorem 16. II we choose to keep the confining external potential fixed. Both 
conventions are equivalent, of course, c.f. the remark at the end of Section 
m but when comparing the exponents of that appear in the proofs of the 
two theorems the different conventions should be born in mind. 

There are two essential components of our proof of Theorem 16.11 The 
first is a proof that the part of the kinetic energy that is associated with the 
interaction v (namely, the second term in (|6.5a|) ') is mostly located in small 
balls surrounding each particle. More precisely, these balls can be taken to 
have radius roughly N~^^^, which is much smaller than the mean-particle 
spacing A^~^/'^. This allows us to conclude that the function of x defined for 
each fixed value of X by 

fx(^) = ^^fo{x,X)>0 (6.7) 

has the property that V^f-^ (x) is almost zero outside the small balls centered 
at points of X. 

The complement of the small balls has a large volume but it can be a weird 
set; it need not even be connected. Therefore, the smallness of V gf^{x) in 
this set does not guarantee that fj^{x) is nearly constant (in x), or even that 
it is continuous. We need f-g{x) to be nearly constant in order to conclude 
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BEC. What saves the day is the knowledge that the total kinetic energy of 
fxi^) (including the balls) is not huge. The result that allows us to combine 
these two pieces of information in order to deduce the almost constancy of 
/j5^(x) is the generalized Poincare inequality in Lemma 16.51 

Using the results of Theorem 16. 31 partial integration and the GP equation 
(i.e., the variational equation for (f*^^ , see Eq. (|5.5j) ) we see that 

lim [ \(p^^{x)\'^\Vgf^{x)\'^dxdX = inNas [ \^^^{x)\^dx . (6.8) 

N~*oo J ^ J 

The following Lemma shows that to leading order all the energy in (|6.8|) is 
concentrated in small balls. 

Lemma 6.4 (Localization of the energy). For fixed X let 



mm |f - Xk\ > iV-V3-5 ^g_g^ 



for some < 6 < 2/9. Then 

lim [dX [ dx\ v^GP (x) 1 2 1 /x (^) I ' = . 



Remark. In the proof of Theorem l4.1l we chose (5 to be 4/51, but the following 
proof shows that one can extend the range of 6 beyond this value. 

Proof. We shall show that 



dX I dx\ip^^{x)\^\Vsf^{x)\' 

X 



k>2 



-SirNa J dX J |(^^^(f)|*|/^(x)|^ 

>-iTTNaj\^^^{x)\^dx-o{l) (6.10) 

as — > CO, which implies the assertion of the Lemma by virtue of 1)6.81) and 
the results of Theorem 16.31 Here, Q% is the complement of The proof 

^^^^^ ^ 

of 1)6.10(1 is actually just a detailed examination of the lower bounds to the 
energy derived in (LSeYlj and |LYlj and described in Sections 2 and 5. We 
use the same methods as there, just describing the differences from the case 
considered here. 
Writing 

fxi^) = llv''^{xk)F{x,X) (6.11) 

k>2 
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and using that F is symmetric in the particle coordinates, we see that 
is equivalent to 



N 

where Qs is the quadratic form 

N , N 



\Q5{F) > -47rNa J |(^GP|4 _ ^(i)^ (6.12) 



QsiF) = [\ViF\^fl\ip'^^{xk)\^dxk 
i=i -^^1 k=i 

N 

+ E vi\x^-x,\)\F\^Yl\^''''{xk)\^dxk 

N 

F|2 J]|^GP(xfc)|2dffc. (6.13) 



i<i<j<N k=l 
N ^ N 



i=l " k=l 

Here fi? denotes the set 

n^i = {{xi,X) G M^^l min|xi - Xk\ < iV"^/^"^}. 

While H6.12|) is not true for all conceivable F's satisfying the normalization 
condition 

N 

/ |F(X,X)|2 J]|^GP(xfc)|2dffc = l, 

k=l 

it is true for an F, such as ours, that has bounded kinetic energy 1)6. 8p . 
Looking at Section 5, we see that Eqs. (|5. 24^ - 1)5. 25|) . H5.44|) - ()5.48|) are sim- 
ilar to (|6.12p . 1)6. 13(1 and almost establish 1)6.12(1 . but there are differences 
which we now explain. 

In our case, the kinetic energy of particle i is restricted to the subset of 
in which min^^j \xi — x^l < N^^^^^^ . However, looking at the proof 
of the lower bound to the ground state energy of a homogeneous Bose gas 
discussed in Section 2, which enters the proof of Theorem 15.11 we see that 
if we choose 6 < 4/51 only this part of the kinetic energy is needed for the 
lower bound, except for some part with a relative magnitude of the order 
e = 0{N~'^°') with a = 1/17. (Here we use the a priori knowledge that 
the kinetic energy is bounded by ()6.8|) . See also the analogous discussion in 
Sectional p. |231) We can even do better and choose some 4/51 < 6 < 2/9, 
if a is chosen small enough. (To be precise, we choose /3 = 1/3 + a and 
7 = 1/3 — 4a in the notation of (|2.56() , and a small enough) . The choice of 
a only affects the magnitude of the error term, however, which is still o(l) 
as — > cxD. 

Proceeding exactly as in Section 5 and taking the differences just men- 
tioned into account we arrive at ()6.12|) . ■ 

In the following, /C C denotes a bounded and connected set that is 
sufficiently nice so that the Poincare-Sobolev inequality (see |LLo) . Thm. 
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8.12) holds on /C. In particular, this is the case if /C satisfies the cone 
property |LLoj (e.g., if /C is a ball or a cube). The next Lemma generalizes 
Lemma 14.21 

Lemma 6.5 (Generalized Poincare inequality). For d>2 let K (ZW^ 
he as explained above, and let h be a bounded function with j]^h = 1. There 
exists a constant C ( depending only on fC and h ) such that for all measurable 
sets O C /C and all f £ H^{IC) with fhdx = 0, the inequality 



\f{x)\^dx<C UjVfix)\^dx+ i^-^Y' l^\Vfix)\'dx\ (6.14) 



IK 

holds. Here \ ■ \ denotes the measure of a set, and Q'^ = IC\Q. 

Proof. By the usual Poincare-Sobolev inequality on fC (see |LLoj . Thm. 
8.12), 

< C\\Vf\\ 

< 2C (||V/||2,,/(,^,)(^^ + ||V/||2,,/(,^,)(^,^) , (6.15) 
if d > 2 and jj^ fh = 0. Applying Holder's inequality 

||V/|L2d/(d+2)(f,) < ||V/||i2(f,)|0|V'^ 

(and the analogue with 0, replaced by Q^), we see that ()6.14() holds with 
C = 2|/C|2/'^(7. ■ 

The important point in Lemma 16.51 is that there is no restriction on 17 
concerning regularity or connectivity. 

Proof of Theorem \6.1\ For some i? > let /C = {f G M?,\x\ < R}, and 
define 

1 

W^^{x)\'^dx Jic 

We shall use Lemma with d = 3, h{x) = \(p^^ {x)\^ / Iv^'^^P, ^ = 
Vt^r\lC and f{x) = f^{x) - {f^)ic (see (jESI and jSH)). Since ip^^ is 
bounded on /C above and below by some positive constants, this Lemma 
also holds (with a different constant C') with dx replaced by \ip^^ {x)\'^dx in 
(|6.14|) . Therefore, 



(fx)>^ = r l,.GP/;^M2^^ / \V>'^''{x)\y-^ix)dx 



dX I dx|v.GP(f)|2[/^(f)-(/^)^]2 
fC 



<C' dX 



\ip'^''{x)\'\V^f^{x)\'dx 
n^nJC 



^25 

E? 



+^ / |^^P(x)P|V.^/^(x)|2<ix 



(6.16) 
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where we used that irj*^ H /C| < (47r/3)A^ . The first integral on the right 
side of H6.16() tends to zero as ^ oo by Lemma 16.41 and the second is 
bounded by (|6.8|) . We conclude, since 



because of the positivity of f-^, that 



1 



liminf — ((^^^|7|99^^) > / |^^^(f)|^dx lim dX dx|*o(^,^)|^ 



Ik. 



n 2 



\ip^^{x)\'^dx 



where the last equality follows from (|5.15p . Since the radius of /C was arbi- 
trary, we conclude that 

1 



lim —{ip 



1, 



implying convergence of -f/N to P'^P in Hilbert-Schmidt norm. Since the 
traces are equal, convergence even holds in trace norm (cf. 'Si', Thm. 2.20), 
and Theorem 16. II is proven. ■ 

We remark that the method presented here also works in the case of a two- 
dimensional Bose gas. The relevant parameter to be kept fixed in the GP 
limit is A^/l ln(a^/)7v)|, all other considerations carry over without essential 
change, using the results in |LSeY2l ILY2j . cf. Sections 3 and 5.2. A minor 
difference concerns the parameter s in Theorem 16.31 which can be shown to 
be always equal to 1 in 2D, i.e., the interaction energy is purely kinetic in the 
GP limit (see |(]Slj ). It should be noted that the existence of BEG in the 
ground state in 2D is not in conflict with its absence at positive temperatures 
|Ho| IM] . In the hard core lattice gas at half filling precisely this phenomenon 
occurs |KLSj . We also point out that our method necessarily fails for the 
one-dimensional Bose gas, where there is presumably no BEG |PiStj . An 
analogue of Lemma 16.41 cannot hold in the ID case since even a hard core 
potential with arbitrarily small range produces an interaction energy that 
is not localized on scales smaller than the mean particle spacing. 



7. The Charged Bose Gas 

The setting now changes abruptly. Instead of particles interacting with a 
short-range potential v{\xi — Xj\) they interact via the Goulomb potential 

^(1^2 Xj I ) — I Xj I 

(in 3 dimensions). There are particles in a large box A of volume as 
before, with p = N/ . 
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To offset the huge Coulomb repulsion (which would drive the particles to 
the walls of the box) we add a uniform negative background of precisely the 
same charge, namely density p. Our Hamiltonian is thus 

N 

H = Y,-^^^^-y{xi)+ <\Xi-Xj\) + C (7.1) 

i=l l<i<ji<Af 

with 

V{x) = p \x — y\^^dy and ^ ~ nP V{x)dx . 

J A 2 

Each particle interacts only with others and not with itself. Thus, despite 
the fact that the Coulomb potential is positive definite, the ground state 
energy Eq can be (and is) negative (just take ^ =const.). This time, large 
p is the 'weakly interacting' regime. 

We know from the work in |LJNj that the thermodynamic limit eo{p) de- 
fined as in (|2.2j) exists. 

Another way in which this problem is different from the previous one is 
that perturbation theory is correct to leading order. If one computes (^', H^) 
with ^ =const, one gets the right first order answer, namely 0. It is the next 
order in l/p that is interesting, and this is entirely due to correlations. In 
1961 Foldy jFj calculated this correlation energy according to the prescrip- 
tion of Bogolubov's 1947 theory. That theory was not exact for the dilute 
Bose gas, as we have seen, even to first order. We are now looking at second 
order, which should be even worse. Nevertheless, there was good physical 
intuition that this calculation should be asymptotically exact. Indeed it is, 
as proved in |LSoj . 

The Bogolubov theory states that the main contribution to the energy 
comes from pairing of particles into momenta k,—k and is the bosonic ana- 
logue of the BCS theory of superconductivity which came a decade later. 
I.e., ^'o is a sum of products of terms of the form exp{ik ■ [xi — Xj)}. 

Foldy's energy, based on Bogolubov's ansatz, has now been proved. His 
calculation essentially implies an upper bound as proved by Dyson in |D2j 
for a slight reformulation of the model. The lower bound is the hard part. 

Theorem 7.1 (Foldy's law). 

p^oo 5r(5/4) \pn J 

This is the first example (in more than 1 dimension) in which Bogolubov's 
pairing theory has been rigorously validated. It has to be emphasized, how- 
ever, that Foldy and Bogolubov rely on the existence of Bose-Einstein con- 
densation. We neither make such a hypothesis nor does our result for the 
energy imply the existence of such condensation. As we said earlier, it is 
sufficient to prove condensation in small boxes of fixed size. 
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Incidentally, the one-dimensional example for which Bogolubov's theory 
is asymptotically exact to the first two orders (high density) is the repulsive 
delta-function Bose gas |LLj . 

To appreciate the —p^^^ nature of (|7.2j) . it is useful to compare it with 
what one would get if the bosons had infinite mass, i.e., the first term in 
H7.1|) is dropped. Then the energy would be proportional to —p^/^ as shown 
in |LNj . Thus, the effect of quantum mechanics is to lower | to |. 

It is supposedly true that there is a critical mass above which the ground 
state should show crystalline ordering (Wigner crystal), but this has never 
been proved and it remains an intriguing open problem, even for the infinite 
mass case. A simple scaling shows that large mass is the same as small 
p, and is thus outside the region where a Bogolubov approximation can be 
expected to hold. 

Another important remark about the — p^/^ law is its relation to the 
_j\r^/^ law for a ^wo-component charged Bose gas. Dyson |D2j proved 
that the ground state energy for such a gas was at least as negative as 
— (const) A^''/^ as A^ ^ oo. Thus, thermodynamic stability (i.e., a linear 
lower bound) fails for this gas. Years later, a lower bound of this —N'^/^ 
form was finally established in |(]LYj . thereby proving that this law is cor- 
rect. The connection of this —N"^/^ law with the jellium —p^/'^ law (for 
which a corresponding lower bound was also given in |CLYj ^ was pointed 
out by Dyson |D2j in the following way. Assuming the correctness of the 
— p^/^ law, one can treat the 2-component gas by treating each component 
as a background for the other. What should the density be? If the gas has a 
radius L and if it has A^ bosons then p = NL~^. However, the extra kinetic 
energy needed to compress the gas to this radius is NL~'^. The total energy 
is then NL~'^ — Np^^^, and minimizing this with respect to L leads to the 
__/\r^/^ law. A proof going in the other direction is in |CLYj . 

A problem somewhat related to bosonic jellium is fermionic jellium. Graf 
and Solovej [(iiSj have proved that the first two terms are what one would 
expect, namely 

eoip) = Ctfp'^' - Cnp'/' + o{p'/^), (7.3) 

where Ctf is the usual Thomas-Fermi constant and Cd is the usual Dirac 
exchange constant. 

As for the dilute Bose gas there are several relevant length scales in the 
problem of the charged Bose gas. For the dilute gas there were three scales. 
This time there are just two. Because of the long range nature of the 
Coulomb problem there is no scale corresponding to the scattering length a. 
One relevant length scale is again the interparticle distance p~^^^. The other 
is the correlation length scale £cor 

~ /3 ' (ignoring the dependence on p). 
The order of the correlation length scale can be understood heuristically as 
follows. Localizing on a scale icov requires kinetic energy of the order of £~^. 
The Coulomb potential from the particles and background on the scale icor 
is (/'^cor)/^cor- Thus the kinetic energy and the Coulomb energy balance 
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when £cor ~ p~^^^- This heuristics is however much too simphfied and hides 
the true complexity of the situation. 

Note that in the high density hmit icor is long compared to the interpar- 
ticle distance. This is analogous to the dilute gas where the scale £c is also 
long compared to the interparticle distance [see p.l2)) ]. There is however 
no real analogy between the scale icor for the charged gas and the scale ic 
for the dilute gas. In particular, whereas eo{p) for the dilute gas is up to 
a constant of the same order as the kinetic energy ~ ju£~'^ we have for the 
charged gas that eo{p) 7^ i^or = p^^'^- The reason for this difference is that 
on average only a small fraction of the particles in the charged gas actually 
correlate. 



7.1. A Short Sketch of the Rigorous Proof. 



Step 1: Localization into small boxes: As mentioned above an im- 
portant step in the rigorous proof is again to divide the big box A into 
smaller boxes of some fixed size i. This time we must require icor ^- We 
again use Neumann boundary conditions on the boundary of each little box. 

In contrast to the dilute gas we can however no longer simply ignore 
the interaction between the small boxes. To overcome this problem we use 
a sliding technique first introduced in |CLYj . To explain this technique 
we introduce the localized interaction w{x,y) = xix)Yi^(x — y)x{y)i where 
Yui{x) = exp(— is a Yukawa potential and x is a smooth ap- 

proximation to the characteristic function of the cube [0,£]^. Consider the 
Hamiltonian 

N 



Hi = '^i-n-f^Aej-p w{xj,y)dy] 
j=l V J J 



+ ^ w{xi,Xj) + lp'^ w{x,y)dxdy, (7.4) 

l<i<j<N 

where —Ai is the Neumann Laplacian for the cube [0, i]^, and 7 is a constant 
depending on x that converges to 1 as x converges to the characteristic 
function of [0,£]^. We consider —Ag as acting in the Hilbert space L^(M^) 
in such a way that it is zero on functions supported away from [0,^]^. Let 
H^^z) be the corresponding Hamiltonian for the translated box z+ [0,^]^. 
If we ignore small errors (that can be fairly easily controlled) due to the 
boundary of A, the original operator H in (|7.H) is bounded below by the 
operator 

rSl^H,{z)dz-^. (7.5) 

The last error term is due to the interaction between the small boxes. This 
requires that the parameter u is chosen appropriately. In fact, it must tend 
to infinity as x converges to the characteristic function of [0, £]^. If we make 
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sure that I » ujp~^/'^ = w^cor we have that uN/l <C Np^^^ and we may 
ignore this error term. 

In the treatment of the dilute gas it was important to understand the 
distribution of particles into the different cells. This is not an issue here. 
In fact the presence of the background will imply that the smallest energy 
is achieved when the number of particles in each cell essentially neutralizes 
the background. More precisely we may simply look for a lower bound on 
the energy independent of the numbers of particles in each little box. 

The operator was defined as an operator on L^(R^^), but we may 
consider its restriction to the invariant subspace with precisely n particles 
in the box [0,£]^. This restriction is equivalent to an operator acting 
on the space L^([0,£]^"). Note that has the form 1)7. 4|) with replaced 
by n, except that now the Neumann Laplacian is considered as acting in 

The problem of finding a lower bound on the ground state energy of H has 
been reduced to finding a lower bound on the operator if" independently 
of n for < n < A^. Without discussing the proof of this point further we 
shall from now on use that it is essentially enough to consider n = pi^ (the 
neutral case). 

Step 2: Reducing to quadratic Hamiltonian: The next step is to 
use the second quantization formalism, which is the one used by Bogolubov 
and which is a very convenient bookkeeping device (but it has to be noted 
that it is no more than a convenient device and it does not introduce any 
new physics or mathematics). In second quantized form we may write the 
Hamiltonian as 

He = P7'^ Yl IpI^^^p + ^ 2Z ^vq,P-i''4^*q^>^^fi ~ "^Op,o/^^g 

+\p"^'^mo^ (7-6) 

where a*, is the creation operator for the eigenfunction u^; of the Neumann 
Laplacian and 

"^pq^flu = jj w{x, y)uf!{x)u^{y)upi{x)up{y) dx dy. (7.7) 

Through a very complicated bootstrapping procedure one now proves that 
to the order of interest here one may ignore several terms in and consider 
instead the Hamiltonian 

p 

pq^O 

We have denoted this operator with a Q since it is quadratic in the operators 
at, with p^O. 
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-Ai<P) > {(p,XeF{-A)Xi(p) , where F{v) = {l-Ct)-^^^jj^, (7.9) 



In order to reduce to this quadratic Hamiltonian it is important to be 
able to control the degree of condensation into the Neumann ground state, 
the constant function tig = More precisely, if we denote by n_|_ = 

a*«p. the operator counting the number of particles not in the con- 
densate, we would like to know that for the minimal energy state, the ex- 
pectation (n+) is small compared to the total particle number n. This, it 
turns out is not too difficult. One needs, however, also a good bound on 
(n^) and this is more difficult. In jLSoj this is not achieved directly through 
a bound on (n^) in the ground state. Rather it is proved that one may 
change the ground state without changing its energy very much, so that it 
only contains values of n+ localized close to This technique, which in 

|LSoj was called "localization of large matrices" , allows one to consider only 
states where (n^) w (^+)^- 

Step 3: Controlling the kinetic energy: The final reduction of the 
Hamiltonian before one can apply the Bogolubov-Foldy method concerns 
the kinetic energy. In fact, if A denotes the Laplacian on we shall use 
the bound (see Lemma 6.1 in |LSoj ) 

'v + {et^y 

for functions <j) orthogonal to constants and where we have assumed that the 
parameter t is chosen such that ||i9'^xl|oo < for all multi-indices a 

with \a\ < 3. Here C > is some universal constant. 

If we introduce the operator bt = a*(xg)ag, where projection of 

the function xi^) exp(iA;-x) onto the orthogonal complement of the constant 
functions, we can write the above inequality as 

J2 Ipf > (2vr)-3n"i I F{\kf)btb^^dk. (7.10) 

p 

Note the extra factor which is due to the fact that b*. contains a factor 

k 

and that we are considering states with total particle number n. 

The reason that we have included a factor in 6*. is in order to be able 

k 

to write the the last sum in Hq also in terms of the 6i. In fact, if we use 
that w{x,y) = xi^)Y^jj{x — y)w{y) we get that 

where 

h{k) = i,r'n~^F{\k\^)(b*h^^ + b*J^^^j 

+ Uk)i-' (btb^^ + b*^^b^^^ + + b^b_^^ . (7.12) 

The last term in (|7.11j) comes from the fact that one must commute the 
factor a^aX in 616- so that it occurs in the normal ordered way aia^ as in 
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Hq. It is not difficult to see that one may ignore the last error term in (|7.11|) 
to the order of interest here. 

Step 4. Applying Bogolubov's method: The last step in the analysis 
is now to use the following simple form of Bogolubov's method. 

Theorem 7.2 (Simple case of Bogolubov's method). 

For constants A> B > and all k we have the operator inequality 

Aibtb^^ + b* ,6 ^) + Bibtb* ^ + b^b 

Proof. We may complete the square 

A{b*^^^ + b*^b ) + Bib*.b*^^ + b^b ^) 

= D{bt + ab ^){bj: + ab* ^) + D{b* + ab^)(b + ab* ^) 
-Do?{[b^^,bV, + [b_^^,b*_J^), 

where 

D{l + a^)=A, 2Da = B. 
We choose the solution a = A/B — A? jB"^ — 1. Hence 

Do? = Ba/2 = ^{A- VA'^-B'^). 

■ 

It is not difficult to see that < ni^. Using the theorem above we 

k k 

see that the ground state energy of i/" up to the errors we have ignored 
is bounded below by -\{2tt)-^ J f{k) - {f{kf - g{kfy/^dk, where g{k) = 
nY^{k) and f{k) = g(k) + f F {\k\'^) . 

Up to the order of interest we may now replace ^~^F{\k\'^) by \k\'^ and 
by the Coulomb potential and thus l^(fe) by 47r|A;|~^. If we now also replace 
n by pl'^ we see that the ground state energy of H"^ to leading order is given 
by 

= 2-1/^-3/4^^3(^/^)1/4 / 1 ^ ^4 _ ^2(2 + ^4)1/2^^^ 

If we finally use that the integral above is 

23/4y^r(3/4) 
5r(5/4) 

and that the leading order of the ground state energy of the original operator 
H according to ()7.5() is I? j^^ times the ground state energy of (again 
ignoring 7) we arrive at Foldy's law (|7.2|1 . 
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